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ESTIMATING  EVAPOTRANSPIRATION 
UNDER  NONHOMOGENEOUS  FIELD  CONDITIONS '  C  0 


By 

Leslie  H.  Parmele  and  Earl  L.  Jacoby,  Jr. 


SUMMARY 


Detailed  micrometeorological  measurements  were 
made  during  selected  days  throughout  1969  by  using 
a  weighing  lysimeter  on  an  adjacent  watershed  that 
had  been  planted  with  corn  in  order  to  test  the 
application  of  commonly  used  meteorological 
methods  of  estimating  evapotranspiration  under  con- 
ditions of  rolling  terrain.  The  experimental  site  was 
selected  on  a  watershed  that  had  a  slope  of  1 3  percent 
ant  that  consisted  of  rolling-surface  topography 
located  in  the  U.S.  Department  of  Agriculture's 
North  Appalachian  Experimental  Watershed, 
Coshocton,  Ohio. 

With  average  hourly  data  being  used,  the  Bowen 
ratio  energy  budget  method  proved  superior  to  the 
stability-adjusted  aerodynamic  methods  in  com- 
paring hourly  lysimeter  flux  rates. 

The  Bowen  method  appears  to  be  capable  of 


accurately  estimating  actual  evapotranspiration 
under  natural  field  conditions  where  fetch  and 
surface  homogeneity  is  less  than  ideal.  Hourly  com- 
putations in  the  combination  equation  for  potential 
evapotranspiration  agreed  with  the  lysimeter  for  bare 
soil  and  for  corn  75  centimeters  high. 

However,  using  the  log  profile  wind  function  in  the 
combination  equation  was  unsatisfactory  for  corn  at 
its  full  height  of  240  centimeters  unless  the  roughness 
length  parameter  (z0)  was  adjusted  as  a  calibration 
term.  An  unusually  low  value  for  z0  of  1 .0  centimeter 
resulted  in  the  best  fit  of  the  combination  equation 
with  lysimeter  flux  rates  for  corn  higher  than  75 
centimeters.  Techniques  were  presented  for  making 
probable-error  analysis  of  evapotranspiration 
estimates  through  meteorological  methods. 


INTRODUCTION 


The  evaporation  from  bare  soils  and  the 
evapotranspiration  from  vegetated  surfaces  have 
been  studied  in  detail  using  meteorological  methods 
capable  of  providing  accurate,  short-period  flux 
rates.  Under  ideal  conditions,  meteorological 
methods  have  been  highly  successful  in  providing 
detailed  information  on  the  vapor  transfer  processes. 


However,  these  methods  have  not  been  tested 
adequately  for  the  nonhomogeneous  surfaces,  soils, 
and  microclimatic  conditions  encountered  in  many 
natural  situations  apart  from  the  ideal  infinite  fetch 
and  short-growing,  well-watered  crops  that 
characterized  most  past  verification  studies. 
Experimental    studies    are    needed    to   test  the 
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application  of  meteorological  theory  against  the 
realities  of  evapotranspiration  from  the  nonideal 
cropped  surfaces  that  are  more  commonly 
encountered  in  agricultural  practices. 

Thus,  the  available  meteorological  methods  should 
be  tested  under  the  realities  of  evapotranspiration 
from  cropped  and  heterogeneous  surfaces.  The 
instrumentation  requirements  and  assumptions  basic 
to  a  given  model  can  be  verified  by  comparing  com- 
puted results  against  an  independent  measure,  pre- 
ferably an  accurate  weighing  lysimeter.  The  purpose 
of  this  publication  is  to  present  such  experimental 
results  under  less-than-ideal  site  conditions. 


Energy  Balance  Method 

One  of  the  most  promising  meteorological 
methods  available  is  an  energy  balance  approach 
known  as  the  Bowen  ratio  method  (2). 3  The  Bowen 
ratio,  defined  as  the  ratio  of  sensible  heat  flux  (H)  to 
latent  heat  flux  (E),  can  be  derived  by  dividing  the 
eddy  transport  equation  for  sensible  heat  flux  by  that 
for  latent  heat  flux. 


B  = 


H  cpPKh  (3T/9Z)p 
E  LeKv(3e/aZ)p 


(1) 


Assuming  Kh  =  K  v,  an  approximation  of  B  can  be 
obtained  by  using  a  finite  difference  form: 


B  =  cpP  (AT)  =  T  (AT) 


L  |  (Ae) 


Ae 


(2) 


where  T  and  e  are  temperature  and  vapor 
pressure  measured  over  the  same  height  interval  just 
above  the  vegetation  or  bare  soil  surface. 

Estimates  of  E  using  the  Bowen  ratio  approach  are 
then  calculated  from: 


R    +  G 

E  =  _J}  


1  + 


(3) 


where  R  n  and  G  are  the  measured  net  radiation  and 
soil  heat  flux  rates.  Equation  3  has  been  used  to 
measure  evaporation  from  homogeneous  surfaces, 
assuming  that  no  horizontal  divergence  of  fluxes 
occurs  between  the  levels  of  measurement  and  the 
evaporating  surface. 


The  advantages  and  disadvantages  of  the  Bowen 
method  are  discussed  in  other  sources  (5, 26).  Tanner 
(27)  reported  close  agreement  between  calculated 
evapotranspiration  rates  and  lysimeter-measured 
values  in  a  humid  climate  with  a  well-watered  crop.  In 
warm,  arid  regions,  the  Bowen  method  has  been  used 
by  Fritschen  (5,  6)  successfully  to  determine  the 
evaporative  flux  for  short  periods.  These  deter- 
minations can  also  be  made  for  longer  period  values. 


Aerodynamic  Approach 

Aerodynamic  approaches  for  estimating  total 
water  vapor  flux  from  an  evaporating  surface  depend 
on  assumptions  that  are  well  documented  in  the 
literature  cited  (23,  26,  28,  32).  Generally,  these 
assumptions  are: 

•  Mean  wind  speed  varies  logarithmically 
with  height. 

•  Eddy  shear  stress  is  constant  with  height. 

•  Eddy  diffusivity  for  water  vapor  (Kv  )  is 
equal  to  that  for  momentum  (Km). 

The  equation  based  on  these  assumptions  most 
commonly  used  in  evaporation  measurement  under 
adiabatic  conditions  is: 


E  =  eL  Pk2(el-e2)  (u2"ul) 


[In  (Z2-D/Zj  -D)] 


(4) 


'Italic  numbers  in  parentheses  refer  to  Literature  Cited,  p.  19. 


Equation  4  permits  the  prediction  of  latent  heat  flux 
rates  from  carefully  measured  profiles  of  vapor 
pressure  (e)  and  horizontal  wind  velocity  (u).  The  dis- 
placement height  (D)  (D  =  z0  +  d)  must  be 
subtracted  from  the  height  value  (Z)  as  measured 
from  ground  surface.  This  adjustment  is  necessary 
because,  as  figure  1  shows,  where  there  is  a  tall,  aero- 
dynamically  rough  vegetation,  the  logarithmic  wind- 
profile  distributions  on  which  equation  4  is  based  are 
displaced  upward  from  the  ground  surface.  The 
parameter  z0  is  a  constant  that  is  used  to  define  the 
roughness  characteristics  of  a  given  surface.  It  is 
determined  by  graphically  extrapolating  the  log- 
arithmic profile  of  measured  velocities  downward  to 
where  wind  speed  theoretically  should  be  zero.  The 
parameter  d  is  the  height  of  the  z0  reference  plane 
displacement. 

The  above  approach  is  valuable  only  under  near- 
adiabatic  conditions.  When  surface  heating  or 
cooling  produces  a  nonadiabatic  temperature  profile, 
the  straight-line  relationship  of  u  versus  1  n  (Z-D)  no 
longer  holds,  especially  under  low-wind  conditions. 
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Figure  1 . — Schematic  representation  of  wind  velocity  distribution 
(U)  with  height  above  ground  (Z)  in  a  plant  community  with 
height  (h).  Wind  velocity  distribution  is  logarithmic  with  height 
above  plants,  and  exponential  with  height  within  the  plant 
community.  Logarithmic  profile  extrapolation  downward  to 
zero  wind  speed  (u  =  o)  is  indicated  by  a  dashed  curve  with  an 
intercept  of  d  +  Z0(I2). 


Ri 


g  (9T/3Z) 
TabsOu/9Z) 


(6a) 


where  g  is  the  constant  of  gravity,  and  where  Tabs  is 
the  absolute  temperature.  Equation  6a  can  be 
approximated  by 


Ri 


(g(Z2-Z1)/Tabs)[(T2-T1)/(u2-u1)2] 


(6b) 


and  is  applicable  to  a  mean  height  (z)  between  z  j  and 
z2,  where  AT  and  Au  are  measured. 

The  coefficient  b  in  equation  5  must  be  determined 
experimentally.  Panofsky,  Blackadar,  and  McVehil 
(16)  recommended  a  value  of  18  after  analysis  of 
several  data  sources.  With  this  b  value,  equation  5  is 
commonly  called  the  "KEYPS"  profile. 

Another  empirical  diabatic  profile  correction 
developed  at  Davis,  Calif.  (20)  is: 


\p  =  (1  ±  50  Ri)  ±  1/4, 


(7) 


where  the  -50  and  -1  /  4  apply  in  the  unstable  case,  and 
the  +  50  and  +1/4  apply  in  the  stable  case. 

Equation  4  can  be  modified  by  introducing  the 
diabatic  correction  function  \Jj  .  Thus,  to  correct 
evaporation  flux  calculations  for  all  conditions  of 
thermal  stability,  the  following  equation  (26)  is  used 


Under  conditions  of  appreciable  temperature 
gradients  because  of  surface  warming  or  "free  con- 
vection," equation  4  will  underestimate  the  flux  rate 
since  the  added  effect  of  buoyant  turbulence  has  not 
been  considered.  Similarly,  temperature  inversions 
over  cold  surfaces  damp  out  turbulent  exchange, 
causing  equation  4  to  overestimate  exchange  rates. 

Several  means  have  been  proposed  to  adjust  for  tne 
influence  of  temperature  gradients  (23).  One 
approach  is  to  incorporate  a  form  of  diabatic 
correlation  function  that  can  then  be  related 
empirically  to  a  stability  parameter. 

The  form  of  the  relationship  normally  used  is: 

^  =  (1-b  Ri)-1/4. 

(5) 

where  Ri  is  the  gradient  Richardson  number,  a  com- 
monly used  stability  parameter.  This  non- 
dimensional  number  (Ri)  is  a  quantitative  ratio  of 
thermal  turbulence  effects,  compared  with  frictional 
turbulence  effects  that  have  the  form  of 


E  =  eLpk    (e  j  -  e? )   (u  2  -  u  j  ) 

[  In  (Z2  -  D)  (Zj  -  D)  ]  2  («//)2  (g) 


Equations  for  Potential  Evapotranspiration 

Potential  evapotranspiration  is  a  micrometeoro- 
logical  concept  that  is  proposed  to  be  the  actual 
evapotranspiration  under  conditions  of  freely  avail- 
able water  and  that  is  intended  to  provide  a  base  line 
for  computing  actual  evapotranspiration  as  a 
fraction  of  potential  evapotranspiration. 

The  combination  method  for  estimating  potential 
evaporation  was  developed  by  Penman  from  both 
energy  balance  and  aerodynamic  transport  equations 
( 17).  Much  attention  has  been  given  to  the  derivation 
and  application  of  the  combination  method  and  its 
various  versions  in  estimating  potential  evaporation. 
Many  sources  are  available  for  more  detailed  infor- 
mation (7,  3,  10,  11,  18,  24,  27,  28,  31). 
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Van  Bavel  (31)  described  the  development, 
limitations,  and  adequacy  of  the  Penman  com- 
bination concept  in  considerable  detail.  More  recent 
formulations  of  the  combination  equation  are  an 
improvement  over  the  original  Penman  versions  in 
that  they  contain  fewer  empirical  constants  or 
functions.  Basically,  the  current  accepted  equation 
(57)  is: 


Ep  =  s       [(Rn  +  G)+'cp     1  (es-ea) 


s  +  7 


s  r. 


in  aerodynamic  transport  characteristics  of  the  layer 
through  which  water  vapor  moves  away  from  the 
evaporating  surface,  assuming  that  vapor  pressure 
(es)  is  at  saturation. 

In  this  form,  no  provision  is  made  for  differences  in 
stomatal  characteristics  making  one  crop  "wetter" 
than  another,  even  when  well  watered.  Tests  of  com- 
bination equations  9  and  1 0  for  open  water,  wet  bare 
soil,  and  well-watered  alfalfa  (a  "wet"  crop)  showed 
excellent  agreement  on  both  an  hourly  and  on  a  daily 
basis  under  a  variety  of  physical  conditions  in 
Phoenix,  Ariz.  (31). 


(9) 


where  e  s  and  ea  are  the  saturated  and  actual  vapor 
pressure  of  the  air  at  height  Za  above  the  given 
surface.  In  contrast  with  the  original  Penman 
formula: 

•  Net  radiation  (Rn  )  is  actually  measured,  instead 
of  relyingon  an  empirical  estimateof  this  parameter. 

•  Correction  is  made  for  soil  heat  flux  (G)  that  had 
been  previously  omited. 

•  r0  is  derived  from  the  log  wind-profile  theory  of 
Thornthwaite  and  Holzman  (29)  and  referenced  from 
the  displacement  height  (D). 


(10) 


•  An  empirical  constant  is  not  included  as  an 
adjustment  for  crop  wetness  or  stomatal  control. 

•  Vapor  pressure  at  the  evaporating  surface  is 
equal  to  the  saturation  vapor  pressure  at  surface 
temperature. 

The  combination  method  for  potential  ET  requires 
these  assumptions: 

1 .  The  eddy  diffusivity  for  momentum  equals  that 


for  heat  and  water  vapor,  K, 


Kv. 


2.  The  air  temperature  at  point  of  measurement 
can  be  substituted  without  undue  error  for  the 
average  of  the  air  and  surface  temperatures  in 
evaluating 

7(S/S  +  7) 

By  referring  to  equation  9,  one  canapply  a  method 
based  on  potential  evapotranspiration  to  any  surface 
with  freely  available  water.  The  only  differences  in 
evaporation  from  a  water  surface,  from  bare  wet  soil, 
and  from  a  well-watered  vegetative  surface  would  be 
due  to  variations  in  surface  radiative  properties  and 


Limitations  of  Meteorological  Methods 

For  routine  field  conditions,  these  meteorological 
methods  could  be  impractical  because  of  measure- 
ment difficulties  and  inaccurate  because  of  the  failure 
of  the  theoretical  model  to  fit  less-than-idea!  con- 
ditions of  a  horizontally  uniform,  two-dimensional 
surface.  Peterson  (19)  discussed  the  problems 
associated  with  flows  encountered  in  heterogeneous 
terrain.  His  predictions  indicate  that  the  standard 
techniques  for  estimating  the  eddy  diffusivities  or 
exchange  coefficients  may  lead  to  unacceptable 
errors  in  estimating  the  vertical  fluxes  of  momentum, 
heat,  and  moisture  over  a  heterogeneous  surface. 
Many  agricultural  and  natural  surfaces  simply  can- 
not be  regarded  as  homogeneous. 

Surface  nonhomogeneity  may  not  only  affect  the 
ratios  of  (Kh  /Km)  and  (Kh  /  Kv  ),  but  also  may  create 
sampling  problems  in  measurement  of  gradients  and 
other  terms  with  the  required  accuracy. 

When  sensors  are  elevated  to  obtain  better  spatial 
samplingand  maximum  possible verticaldifferences, 
the  profile  measurements  of  wind,  temperature,  and 
vapor  concentration  are  affected  more  and  more  by 
surfaces  further  upwind  and  less  by  the  surface 
immediately  below  the  sensor.  At  some  point  the 
measured  profile  gradients  are  not  representative  of 
processes  that  are  occurring  on  the  reference  surface. 
Thus,  a  buffer  area  of  the  same  type  as  the  reference 
surface  is  required  to  reduce  the  possibility  of 
divergence  errors  in  measuring  the  vertical  fluxes. 
This  buffer  area  is  referred  to  as  "fetch"  and  is  the 
minimum  upwind  distance  to  a  significantly  different 
surface  from  the  point  of  measurement. 
Nonhomogeniety  influences  fetch  requirements  in 
that  the  required  height-to-fetch  ratios  should  be  1  to 
50  for  minor  discontinuities  and  at  least  1  to  100  or  1 
to  200  if  the  discontinuity  is  severe  (26).  These  "ideal" 
fetch  requirements  are  difficult  if  not  impossible  to 
satisfy  under  many  field  situations. 
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EXPERIMENTAL  PROCEDURE 


The  experimental  site  was  developed  on  water- 
shed 109  at  the  U.S.  Department  of  Agriculture, 
North  Appalachian  Experimental  Watershed. 
Coshocton,  Ohio.  The  watershed  (fig.  2)  lies  just 
below  the  crest  of  a  major  ridge  that  extends  through 
the  area  and  that  has  an  average  slope  of  1 3  percent  in 
an  eastward  direction.  Agricultural  fields  border  the 
north,  south,  and  west  sides  of  the  watershed;  gravel 
roads  traverse  the  north  and  east  sides.  The  rolling 
terrain  is  far  from  an  ideal  homogeneous  surface  as 
required  for  the  classical  type  of  micro- 
meteorological  studies.  At  the  center  of  the  water- 
shed a  micrometeorological  sampling  site  for  temper- 
ature, humidity,  and  wind  profiles  was  established. 
Most  energy  budget  data  and  further  profile 
measurements  of  humidity,  temperature,  and  wind 
were  obtained  at  lysimeter  Y102C.  located  200  feet 
north  of  the  central  watershed  sampling  site.  The 
instrument  trailer  that  housed  the  recording  equip- 
ment was  parked  on  the  boundary  of  the  field 
approximately  100  feet  northwest  of  the  lysimeter 
site. 

The  entire  watershed,  including  lysimeter.  was 
planted  in  corn  on  May  1,  1969.  in  42-inch  rows 
generally  oriented  in  a  southwest-northwestdirection 
parallel  to  the  contour  of  the  field.  The  prevailing 
winds  in  the  area  are  from  the  southwest  and  provide 
the  best  fetch  conditions  for  the  two  sites. 

Complete  micrometeorological  measurements 
were  made  at  the  lysimeter  site  on  May  2  and  22. 
June  26.  July  24,  August  27,  and  October  14and  16, 
1969.  Micrometeorological  measurements  were  not 
made  at  the  midwatershed  site  during  the  first  two 
test  periods.  The  midwatershed  profile  data  on 
June  26  were  inaccurate  because  of  instrumental 
difficulties,  so  that  data  at  the  midwatershed  site  were 
available  only  on  July  24  and  afterward.  A  summary 
of  conditions  for  each  data  is  given  in  table  1 . 

The  corn  was  harvested  on  October  1 5  and  the  soil 
disked  to  prepare  for  winter  wheat  planting.  Thus, 
the  data  on  October  16  were  obtained  on  bare  soil 
partly  covered  with  litter  from  cornstalks. 

A  complete  description  of  the  instrumentation  and 
data-reduction  techniques  can  be  found  in 
appendix  B.  A  listing  of  the  average  hourly  micro- 
meteorological data  is  listed  in  appendix  C. 


Variation  of  D  and  Z  0 

The  determinations  of  D  and  Z  0  that  were  used  in 
equations  4  and  10  are  two  of  the  most  demanding 
measurements  in  the  profile  methods,  because  each 
depends  on  accurate  measurement  of  the  second- 
order  properties  of  the  profile.  D  and  Z0  must  be 
determined  under  adiabatic  situations  because  lack 
of  stability  in  these  paramenters  may  introduce  errors 
(15).  Since  D  and  Z0  of  cropped  surface  may  change 
several  centimeters  with  varying  wind  speed  (12,  30), 
errors  may  also  arise  from  assuming  that  D  and  Z0 
are  constant  for  a  given  surface. 

Mukammal  (15)  has  indicated  that  the 
determination  of  D  is  very  sensitive  to  small  errors  in 
wind  speed  measurement.  He  observed  that  under  a 
given  set  of  conditions,  an  increase  of  only  1.0 
centimeter  per  second  in  the  average  speed  of  the 
middle  level  of  a  three-level  wind  system  increased  D 
by  50  percent,  from  40  centimeters  to  60. 

Error  Analysis 

The  meteorological  estimates  of  E  were  basically 
indirect  measurements:  that  is,  the  unknown  (E)  was 
determined  by  using  an  equation  that  related  it  to 
some  other  measurable  quantities  such  as  radiation, 
temperature,  wind  speed,  and  vapor  pressure.  Each 
measured  variable  had  a  certain  standard  deviation 
of  the  mean  or  probable  range  of  error.  Thus,  each 
computed  value  of  E  had  some  uncertainty  associated 
with  it.  The  magnitude  of  the  uncertainty  depended 
both  on  the  magnitude  of  the  errors  in  each  of  the 
measured  variables  and  on  the  nature  of  the  equation 
itself. 

Because  of  the  problems  in  evaluating  D  and  Z 
accurately,  the  probable  maximum  error  for  each  was 
assumed  to  be  ±  25  percent  of  the  absolute  value  for 
use  in  error  analysis.  The  magnitude  of  these  possible 
errors  appeared  reasonable  in  light  of  the  preceding 
discussion. 

A  probable  uncertainty  or  error  in  each  calculated 
hourly  value  of  E  was  evaluated  by  using  procedures 
outlined  in  appendix  A.  All  estimates  of  possible 
error  were  based  on  probable  measurement  errors  for 
each  variable  (table  2).  The  estimated  range  or 
probable  uncertainty  in  a  particular  E  estimate 
assumed  that  all  sensors  and  recording  systems  were 
operating  within  their  accepted  limits  and  that  no 
error  was  introduced  in  the  data-reduction  process. 


©  -TEMPERATURE  &  HUMIDITY 
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Figure  2.  —  Map  of  field  experimental  site  at  Coshocton,  Ohio,  showing  location  of  instrument  trailer,  lysimeter  YI02.  and 
micrometeorological  sensors. 
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Table  I.— A  summary  of  conditions  and  measurements  on  days  of  experimental  observations  at 
lysimeter  Y102C  and  watershed  109,  Coshocton,  Ohio 


Soil 

Accumu- 

Accumu- 

Cover 

Av.  u 

General 

Height 

Length 

Date 

Time 

condi- 

lated E 

lated 

and 

(cm/s) 

weather 

Measure- 

of Ae 

of 

(EST) 

tions 

(lysi- 

precipi- 

crop 

direction 

contions 

ments 

and  AT 

fetch 

meter) 

tin 

height  (h) 

Cm 

Cm 

Cm 

Cm 

M 

5/2/69 

04-22 

Dry  and 

0 

0 

Bare 

250  (SW) 

Clear  and 

Lysimeter 

10  and  100 

150 

loose 

soil. 

warm. 

Y102C,Rn, 

D  =  0 

on 

surface. 


G,  Ae, 
AT,  Au. 


=  0.1 


5/22/69    04-22  Surface 
moist. 


5.1 


6.1 


General 
bare  soil 
corn 
h  =  5. 


150  (N) 


Overcast 
and  cold 


Lysimeter 
Y102C,Rr 


0.05    rain     G,  Ae, 
@  1200  EST.  AT,  Au. 
EST. 


10  and  100 
D  =  0 

z„  =  0 


15 


6/26/69    04-22    Wet  soil  15.2 
saturated 
to  field 
capacity. 


20.9 


7/24/69    04-22  Wet 


10/16/69  07-19 


Dry, 

surface 

moist. 


26.7 


8/27/69    05-23    Surface  42.0 
dry,  sub- 
surface 
moist. 


10/14/69    05-20  Dry 


Corn 
almost 
complete 
cover 
h  =  75. 


39.4  Corn 

h  =  240. 


52.0  Corn 

h  =  240. 


200  (SW) 


57.5 


57.7 


Corn 
ready 
to 

harvest 
h  =  130. 


Bare 
soil 
with 
litter. 


100  (SW) 


100  (N-S) 


150  (NW) 


500  (SW) 


Clear  and 
warm. 


Clear  with 
haze,  rain 
@  1500 
EST 

Clear, 
light, 
variable 
wind. 


Overcast 
and  cool. 


Partly 
cloudy 
and  cool. 


Lysimeter 
Y102C.Rr 
G,  Ae, 
AT,  Au. 


Lysimeter 
Y102C,Rn, 
G,  Ae, 
AT,  Au. 

Lysimeter 
Y102C,Rn, 
G,  AE, 
AT,  Au 
Watershed 
109  Ae,  AT. 

Lysimeter 
Y102C,  Rn, 
G,  Ae, 
AT,  -Hi 
Watershed 
109.  Rn 
Ae,  AT, 
Au. 

Lysimeter 
Y102C,  Rn, 
G,  Ae, 
AT,  Au 
watershed 
109,  Rn, 
Ae,  AT, 
Au. 


150  and  250 
D  =  0 
zQ  =  2.0 


250  and  350 
D  =  150 
z=  15.0 


250  and  350 
D  =  150 
zQ  =  20.0 


150  and  250 
D  =  50 


150 


150 


15  to  150 
50  to  75 


15  to 
100 


=  8.0 


25  and  100 
D  =  o 
=  0.1 


150 


Table  2.— Recording  and  sensor  accuracy  of  micrometeorological  variables  for  use  in  probable-error  analysis 


V  til  It!  UlC 

Accuracy  of 

Probable 

m  t^o  c  1 1      m  t 

[llCaaUiCIIldl  I 

Rn 
n 

±0.01  langley/min 

±0.02  langley/min 

±0.02  langley/min 

G 

±  .01  langley/min 

±  .02  langley/min 

±  .02  langley/min 

±  .35°  C 

±  .10°  C 

±  .40°  C 

AT 

±  .02°  C 

±  .02°  C 

±  .03°  C 

AT(wet  bulb) 

±  .02°  C 

±  .02°  C 

±  .03°  C 

es,ea 

±  .60  mbar 

Ae 

*  .10  mbar 

(e.-e.) 
s  a 

*  .80  mbar 

u 

±4.5  cm/s 

--- 

±4.5  cm/s 

A  1 1 

~  f\  A  f  m  /  c 

zo 

Assume  25%  of  ZQ 

D 

Assume  25%  of  D 

EXPERIMENTAL  RESULTS 


Bowen  Ratio  and  Aerodynamic  Methods 

A  comparison  of  hourly  estimates  E  using  the 
Bowen  ratio  and  aerodynamic  equations  with  the 
lysimeter  value  for  M  ay  2  is  shown  in  figures  3  and  4, 
and  for  June  26  in  figures  5  and  6,  respectively.  Both 
days  were  characterized  by  clear  skies.  Wind  was 
from  the  south-southwest  for  most  of  the  day,  with 
speeds  ranging  from  100  to  300  centimeters  per 
second.  This  direction  provided  the  most  ideal  fetch 
conditions  at  the  lysimeter  test  area. 

On  May  2,  measurements  were  made  over  bare, 
dry  soil.  As  a  result  of  the  solar  heating  of  the  bare 
soil,  measurements  indicated  strong  temperature  and 
humidity  gradients,  and  the  soil  surface  temperature 
reached  40°  C  during  midday.  Bowen  ratios  (B) 
ranged  from  1.2  to  1.9  during  daylight  hours, 
indicating  that  most  of  the  radiant  energy  was  being 
transferred  into  sensible  heat. 

The  large  negative  Richardson  numbers  that 
prevailed  during  the  measurements  on  May  2 
resulted  in  a  large  diabatic  correction  to  the 
aerodynamic  equation  (8).  Evaporation  computed 
from  the  simple  uncorrected  logarithmic  equation  (4) 
averaged  50  percent  less  than  those  obtained  with 
equation  8.  This  required  the  largest  diabatic 
correction  on  any  of  the  7  test  days.  The  closeness  of 
the  Bowen  ratio  estimate  and  the  aerodynamic 
estimates  using  "KEYPS"  correction  was 
encouraging.  There  appeared  to  be  a  tendency  for  the 


aerodynamic  equation  with  the  Davis  adjustment  to 
overcompute  E  on  this  day. 

The  range  of  probable  error  on  May  2  was  very 
small  for  both  aerodynamic  and  Bowen  ratio 
estimates,  ranging  from  ±  0.01  to  ±  0.02  langley  per 
minute.  The  probable  error  for  the  Bowen  ratio 
estimate  was  shown  to  increase  considerably  at  1800 
hours  to  a  range  of  ±  0.07  langley  per  minute. 

The  run  of  June  26  was  made  over  corn  75 
centimeters  high  following  a  week  of  heavy  rains  that 
had  brought  the  soil  to  field  capacity.  This  day  was  a 
marked  contrast  to  conditions  on  May  2,  in  that  the 
temperature  profiles  were  very  close  to  adiabatic 
during  most  daylight  hours.  Practically  all  the  net 
radiation  went  into  latent  heat  flux,  with  Bowen  ratio 
values  (B)  ranging  from  -0.13  to  0.20  during  day- 
light hours.  Because  of  the  low  Richardson  numbers, 
(-0.01  to  +0.01 ),  only  a  slight  diabatic  correction  (i/>  2) 
was  necessary  for  the  aerodynamicequation.  ranging 
from  0.90  to  1.0,  maximum. 

Evaporation  flux  rates  computed  by  aerodynamic 
and  Bowen  ratio  methods  on  June  26at  the  lysimeter 
site  were  compared  with  the  measured  lysimeter 
values  (figs.  5  and  6).  The  Bowen  ratio  method 
agreed  closely  with  the  lysimeter  data,  while  the  aero- 
dynamic method  provided  slightly  lower  values 
during  the  hours  1  100  to  1600. 

Except  for  a  couple  of  hours,  the  Bowen  ratio 
estimate  was  generally  within  the  probable  error 
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Figure  3.  —  Hourly  evaporative  flux  from  bare  soil  obtained  on  May  2.  1 969.  from  weighing  lysimeter  Y  102Cat  Coshocton.  Ohio, 
and  calculated  with  Bowen  ratio  method  showing  range  of  probable  error  in  estimating  LE. 


LYSIMETER  SITE 
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-•-LYSIMETER 
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I    SHOWING  PROBABLE  RANGE  OF  ERROR 

-A— AERODYNAMIC  WITH  "DAVIS"  ADJUSTMENT 


10  12         14  16  18 

TIME  OF  PERIOD  (hrs.  E.S.T.) 


20 


24 


Figure  4.  —  Hourly  evaporative  flux  from  bare  soil  obtained  on  May  2.  1969.  from  weighing  lysimeter  Y  1 02C  at  Coshocton.  Ohio, 
and  calculated  with  adjusted  aerodynamic  methods  showing  range  of  probable  error  in  estimating  LE. 
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Figure  5.  -Hourly  evaporative  flux  from  corn  (75  cm)  obtained  on  June  26.  1 969.  from  weighing  lysimeter  Y I02C  at  Coshocton. 
Ohio,  and  calculated  with  Bowen  ratio  method  showing  range  of  probable  error  in  estimating  LE. 
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LYSIMETER  SITE 
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Figure  6.  —  Hourly  evaporative  flux  from  corn  (75  cm)  obtained  on  June  26,  1 969.  from  weighing  lysimeter  Y  102C  at  Coshocton.  Ohio, 
and  calculated  with  adjusted  aerodynamic  methods  showing  range  ot  probable  error  in  estimating  LE. 


range  of  0.03  langley  per  minute  from  the  lysimeter 
value.  The  range  of  probable  error  for  the  aero- 
dynamic estimate  was  very  large  from  1 100  to  1600 
hours,  ranging  from  ±0.10  to  ±0.15  langley  per 
minute  on  an  hourly  basis.  Reliability  in  the  aero- 
dynamic estimate  of  E  was  greater  during  the  early 
morning  and  late  evening  hours. 

D  was  assumed  to  be  accurate  within  ±  25  percent 
of  the  "true"  value  for  use  in  the  error  analyses.  On 
June  26th,  with  h  =  75  centimeters  and  D  -  50 
centimeters,  this  probable  error  in  D  would  be  ±  12 
centimeters.  If  it  were  possible  to  estimate  D  without 
error,  the  range  of  probable  uncertainty  in  the 
aerodynamic  estimate  of  E  would  be  reduced  to  0 .06 
langley  per  minute  from  1  100  to  1600  hours. 

Thus,  the  large  uncertainty  in  the  aerodynamic 
estimate  can  be  traced  partly  to  the  uncertainty  in  the 
estimate  of  D.  Confidence  could  be  improved  by 


increasing  the  sampling  heights  so  that  D  is  small, 
relative  to  Z2  and  Z  1 .  However,  if  this  is  done,  the 
profile  measurements  may  not  be  representative  of 
the  processes  occurring  on  the  given  surface. 

The  comparisons  described  above  are  for  the  2  test 
days  considered  most  ideal.  Conditions  for  the  other 
days  were  often  far  from  ideal.  A  summary  of  total 
daylight  vapor  flux  (E)  by  the  various  methods  for  all 
periods  is  given  in  table  3  and  is  listed  in  appendix  D 
along  with  the  lysimeter  values  for  the  hourly  periods. 

The  Bowen  ratio  estimate  was  within  1 0  percent  of 
the  total  lysimeter  E  measured  on  five  of  the  seven 
cases  and  within  20  percent  on  the  last  2  test  days.  The 
difference  for  the  total  E  on  these  days  is  well  within 
the  accuracy  of  the  meteorological  and  lysimeter 
measurements. 

The  results  obtained  with  the  aerodynamic 
estimates  were  more  variable.  Only  in  a  few  cases  did 
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Table  3.— A  comparison  of  observed  lysimeter  E  ,  with  E  computed  hourly  by  Bowen  ratio  and 
aerodynamic  methods  with  daylight  total  added 


Date 

I  line 
(EST) 

E 

Lysimeter 
Y102C 

Bowen  ratio 

Aerodynamic  equation 

Fetch 
height 

Equation  3 

Probable 
error 

Davis  adj. 
equation 
8  and  5 

Keyps  adj. 
equation 
8  and  7 

Probable 
error 
Keyps 

Mm 

Mm 

Mm 

Mm 

Mm 

Mm 

Lysimeter  site 

5/2/69 

7-18 

1.42 

2 1.40 

±0.08 

2.42 

2 1.57 

--- 

150 

5/22/69 

12-20 

1.23 

2 1.33 

±  .11 

2 1.27 

2 1.29 

15 

6/26/69 

5-22 

7.06 

2  6.74 

±  .12 

26.35 

35.90 

+0.37 

50 

7/24/69 

6-16 

3.81 

2  3.46 

±  .12 

3  3.31 

2.30 

±  .56 

45 

8/27/69 

6-23 

4.60 

24.38 

±  .16 

2.78 

l  .y\j 

±  .17 

5-45 

10/14/69 

8-18 

.92 

.60 

±  .08 

.49 

.38 

±  .04 

6 

10/16/69 

8-17 

1.10 

3  .90 

±  .16 

2 1.22 

2 1.20 

±  .01 

150 

Watershed  site 

7/24/69 

6-16 

3.81 

33.20 

±  .11 

5.76 

2  3.78 

±1.08 

20 

8/27/69 

6-23 

4.60 

34.10 

±  .21 

3.05 

2.10 

±  .37 

15-20 

10/14/69 

8-18 

.92 

2  .94 

±  .17 

3 1.30 

3 1.10 

±  .08 

40 

10/16/69 

8-17 

1.10 

.46 

±  .5 

.08 

.07 

±  — 

75 

1  Ratio  of  minimum  upwind  distance  from  sampling  site  of  surface  being  measured  to  height  of  highest  sampling  point  in  profile 
measurements  used  in  computations. 

2Computed  E  within  ±10  percent  of  lysimeter  E. 
3  Computed  E  within  ±20  percent  of  lysimeter  E. 


the  computed  flux  rates  closely  agree  with  the 
lysimeter  or  Bowen  ratio  estimates  of  E  after  crop 
height  had  exceeded  75  centimeters.  Between  May  2 
and  July  24,  the  probable  range  of  uncertainty 
increased  from  ±  0.07  to  ±  0.56  langley  per  minute,  as 
D  increased  because  of  crop  growth  from  0  to  1  50 
centimeters.  The  total  flux  estimate  on  July  24  would 
have  had  an  uncertainty  of  ±  0.24  langley  per  minute 
if  it  could  be  assumed  that  "D"  had  been  measured 
without  error.  Thus,  a  major  source  of  error  in  the 
aerodynamic  method  appears  to  be  the  uncertainty  in 
the  measurement  or  approximation  of  the  displace- 
ment height  D. 

A  large  variation  occurred  in  the  aerodynamicand 
Bowen  ratio  estimates  at  the  lysimeter  site  on 
August  27  and  again  on  October  1 4.  This  failure  can 
be  related  to  poor  fetch  conditions;  light,  variable 
winds  arose  on  both  days.  On  October  14,  the  water- 
shed site  was  better  suited  in  regard  to  wind  direction 
and  fetch.  The  results  at  the  watershed  site  showed 
that  all  three  estimates  agreed  on  this  day.  In 
contrast,  on  October  1 6  and  under  heavy  winds,  the 
watershed  site  failed  to  provide  any  meaningful 
estimates,  whereas  the  lysimeter  site  with  a  longer 
fetch  closely  agreed  on  all  three  estimates.  Generally, 
the  Bowen  ratio  estimates  of  E  provided  reasonable 


accuracy  when  the  fetch  /  height  ratio  was  in  the  range 
of  20  to  1  through  50  to  1,  except  for  October  16 
when,  for  bare  soil,  a  fetch  of  75  to  1  was  inadequate. 
A  fetch/ height  ratio  of  greater  than  50  to  1  appears 
necessary  for  successful  application  of  aerodynamic 
techniques. 

Combination  Formula 

Four  days— May  22,  June  26,  July  24,  and 
August  27 — were  considered  to  have  met  the  require- 
ments of  nonlimiting  water  supply  necessary  to 
satisfy  the  definition  of  potential  evapotranspiration. 
Latent  heat  flux  was  computed  hourly  for  each  day 
from  equation  9  using  Z0  as  determined  from  wind- 
profile  analysis.  The  results  for  two  of  these  days  are 
plotted  in  figures  7  and  8.  On  June  26,  the  estimated 
value  of  E  p  was  close  to  E,  but  about  10  percent 
greater  than  the  amount  as  measured  by  the 
lysimeter.  The  range  of  probable  error  in  the  Ep  was 
generally  ±0.04  to  ±0.08  langley  per  minute, 
assuming  an  uncertainty  in  Z0  of  0.5  centimeter, 
with  Z0  equal  to  2.0  centimeters.  Thus  considering 
the  range  of  uncertainty,  hourly  Ep  estimates  agreed 
with  the  lysimeter  except,  perhaps,  from  1400 
to  1700. 
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-Hourly  evaporative  flux  from  corn  (75  cm)  obtained  on  June  26.  1 969.  from  weighing  lysimeter  Y  1 02C  at  Coshocton, 
and  calculated  potential  ET  by  combination  equation  showing  range  of  probable  error  in  estimating  II 

LYSIMETER  SITE 
JUNE  26,  1969 
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Figure  8.  —Hourly  evaporative  flux  from  corn  (240  cm)  obtained  on  August  27. 1 969.  from  weighing  lysimeter  Y  1 02C  at  Coshoc- 
ton, Ohio,  and  calculated  potential  ET  by  combination  equation  showing  range  of  probable  error  in  estimating  I.E. 
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In  contrast,  the  results  for  August  27  revealed  an 
overestimation  of  E  by  more  than  1 00  percent  of  the 
lysimeter  value.  Ep  also  exceeded  net  radiation  by  1 0 
to  15  percent  from  1100  to  1600.  The  range  of 
probable  error  (±0.15)  was  very  large  compared  with 
the  June  26  range  in  probable  error  (  ±0.06). 

A  summary  of  totalized  daylight  hourly  values  for 
each  of  the  4  days  is  given  in  table  4,  column  2. 
Equation  9,  using  measured  Z0  values,  failed  to 
predict  E  accurately  once  the  corn  had  reached  full 
height  at  240  centimeters.  On  May  22,  agreement 
between  computed  Ep  and  lysimeter  flux  rates  for 
bare  soil  was  excellent  despite  poor  fetch  conditions. 

The  value  Z0  can  be  estimated  as  a  function  of  crop 
height  (26)  by  the  equation: 

ZD  =  0.058  h  119  (n) 

The  value  Z  0  could  have  been  estimated  on  the  basis 
of  equation  1 1  rather  than  on  the  basis  of  a 
transposed  value  from  the  Ithica,  N.Y.,  experiment 
previously  mentioned  or  wind-profile  analysis.  Then, 
the  values  of  Z  0  would  be  10  and  40  centimeters  for 
crop  heights  of  75  and  240  centimeters,  respectively. 
Used  in  equation  9  these  Z  0  values  would  greatly 
increase  the  difference  between  computed  E  p  and 
observed  E.  Thus,  it  was  necessary  to  substitute  the 
value  from  the  Ithaca,  N.Y.,  experiment  because  of 
the  nonhomogeneous  terrain  at  the  Coshocton,  Ohio, 
site. 

In  reference  to  the  totalized  daylight  values  given  in 
table  4,  the  probable  error  in  the  period  total  for 
individual  days  increased  considerably  as  Z0  rose 


from  0. 1  to  20  centimeters,  assuming  an  uncertainty 
of  ±25  percent  in  the  Z0  estimate.  An  absolute 
measure  of  Z  0  would  reduce  the  range  of  probability 
error  to  ±0.13  millimeter  per  day  or  less  for  all  4  days. 
This  uncertainty  in  the  E  p  estimate  could  account  for 
some  of  the  difference;  however,  it  is  doubtless  not 
the  only  reason  for  the  overcomputation  of  E  on 
July  24  and  August  27. 

The  log  law  for  neutral  profiles  from  which 
equation  10  was  derived  is  recognized  to  be  in  error 
for  diabatic  situations.  Tanner  (26)  suggests  that  the 
wind  function  (10)  includes  a  diabatic  convection 
factor  similar  to  that  used  in  the  aerodynamic 
equation  correction.  However,  if  such  a  correction 
factor  were  applied  to  the  3  days  in  question,  it  would 
increase  the  error  between  estimated  and  observed 
lysimeter  flux  rates.  Accounting  for  this  disagree- 
ment would  be  easy  by  assuming  that  the  surface  was 
not  a  potential  evaporation  surface  or  that  the 
measurements  of  the  meteorological  variables  were 
not  representative.  The  success  of  the  Bowen  ratio 
estimates  invalidates  this  later  argument. 

If  the  surface  is  not  evaporating  at  the  potential 
rate,  then  computing  the  surface  diffusion  resistance 
(rs),  which  is  a  means  of  defining  the  degree  of 
unsaturation  of  the  crop  and  soil  system,  should  be 
possible.  When  the  energy  balance  of  the  crop  and  the 
aerodynamic  resistance  term  are  known,  rs  can  be 
calculated  as  a  residual  from  the  ratio  of  potential  to 
actual  evaporation  (25). 


Table  4.—  Observed  lysimeter  E  during  periods  of  adequate  moisture  supply,  compared  with  potential 
evaporation  (Ep),  computed  hourly  by  combination  equation  with  total  for  period 


E 

Equations 

Probable 

Equations 

Date 

Time 

Lysimeter 

9  and  10 

error 

9  and  10 

Penman 

Measured 

(EST) 

Y102C 

measured  zQ 

zQ  ±25% 

assume 

version 

zo 

zQ  =  1 .0  cm 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

Mm 

Mm 

Mm 

Mm 

Mm 

Cm 

5/22/69 

1200  to  2000 

1.23 

1  1.06 

±0.07 

'  1.17 

0.1 

6/26/69 

500  to  2200 

7.06 

2  7.85 

±  .20 

'  7.14 

'6.78 

2 

7/24/69 

600  to  1600 

3.81 

4.55 

±  .28 

'  3.51 

1  3.60 

15 

8/27/69 

600  to  2300 

4.60 

8.14 

±  .47 

'4.87 

'5.03 

20 

1  Computed  Ep  within  ±10  percent  of  lysimeter  E. 
2 Computed  E   within  ±20  percent  of  lysimeter  E. 
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Potential  evaporation  (Ep  )  can  be  calculated  from 
equation  9.  When  the  measured  Z  0  values  were  used, 
the  calculation  of  rs  provided  hourly  values  from  0  to 
0.8  second  per  centimeter  for  June  26,  and  from  0.5 
to  2.5  seconds  per  centimeter  for  August  27  during 
daylight  periods.  These  estimates  for  rs  appear 
reasonable,  but  there  are  no  independent  measures  of 
rs  or  soil  water  potential  for  comparison.  Because  of 
the  high  degree  of  recharge  during  this  period 
(table  2).  precipitation  -  measured  evaporation), 
probably  no  soil  moisture  deficiency  existed.  Thus, 
water  was  not  limiting,  and  potential  conditions 
should  be  as  valid  on  July  24  and  August  27  as  on 
June  26. 

Jensen  (9)  and  Rosenberg  (21)  both  report  failure 
of  the  combination  equation  9  when  using  the 
aerodynamic  term  (10)  proposed  by  Van  Bavel  (J/). 
Jensen  indicates  that  under  arid  conditions  it  is 
necessary  to  calibrate  ra  by  varying  the  Z0  value 
until  the  E  p  estimate  of  equation  9  had  agreed  with  a 
lysimeter  measure.  Sellers  (22)  has  indicated  that 
unreasonably  high  E  estimates  have  been  obtained  by 
using  simple  logarithmic  profile  methods  over  tall 
crops.  Van  Bavel's  verification  of  the  modified  com- 
bination equation  (31)  was  limited  to  bare  soil  and 
alfalfa  25  to  30  centimeters  high.  The  authors  too 


could  find  no  reported  applications  of  equation  9  to  a 
tall  crop  such  as  corn,  where  an  independent  measure 
of  E  was  available  for  verification. 

Potential  evaporation  was  recalculated  for  the  3 
days  with  crop  cover  byequation9,usinganassumed 
Z0  of  1.0  centimeter  in  the  wind-function  term  ( 10). 
Improvement  in  the  period  totals  for  each  of  the  3 
days  can  be  observed  in  table  4,  column  4.  While  the 
change  in  Z0  from  2  centimeters  to  1  for  June  26 
resulted  in  a  slight  improvement,  the  reduction  in  z.  o 
from  20  centimeters  to  1  on  August  26  resulted  in  an 
Ep  estimate  that  is  within  5  percent  of  the  lysimeter 
value,  as  compared  with  a  previous  overestimate  of 
100  percent.  These  results  agree  favorably  with  those 
of  Jensen  (9),  who  reported  that  for  a  variety  of  crops 
in  southern  Idaho.  Z0  in  equation  10  should  be 
approximately  0.75  centimeter  where  crops  were 
higher  than  50  centimeters. 

For  comparison,  the  combination  equation  9 
using  Penman's  original  aerodynamic  term  (17)  was 
computed  hourly  and  summed  to  give  a  total  for  each 
day.  The  results  are  shown  in  table  4  column  5 .  They 
indicate  that  the  wind  function  as  originally  used  by 
Penman  produces  essentially  the  same  results  as  the 
log  profile  version  (10)  using  a  Z0  value  of  1 
centimeter. 


DISCUSSION 


Most  of  the  meaningful  verifications  of  the  Bowen 
ratio  energy  budget  technique  have  been  reported  for 
low-growing  irrigated  crops  with  generally  flat 
topography  and  uniform  fetch  conditions.  The 
success  of  the  hourly  values  of  E  computed  by 
equation  3  for  the  lysimeter  at  Coshocton.  Ohio, 
tends  to  verify  the  application  of  the  approach  to 
estimating  evapotranspiration  under  natural  field 
conditions  where  fetch  and  surface  homogeneity  are 
less  than  ideal.  The  method  is  capable  of  predicting 
actual  evaporation  from  a  bare,  dry  soil  with  limiting 
moisture,  as  well  as  from  a  tall  crop  such  as  corn 
(height  =  240  centimeters)  under  a  wide  range  of 
growth  and  moisture  conditions. 

Overcast  skys  (low  direct  solar  radiation)  or 
variable  weather  conditions  do  not  appear  to  limit  the 
application  of  this  technique  except,  perhaps,  during 
periods  of  rain.  The  Bowen  ratio  method  fails  during 
early  morning  and  late  evening  as  a  result  of 
increasing  relative  errors  in  measuring  temperature 
and  humidity  profiles  during  nondaylight  hours. 
However,  since  E  rates  are  quite  small  and  are  often 
zero  at  such  times,  this  failure  is  relatively 
unimportant. 


Use  of  average  hourly  values  of  AT  and  Ae,  as 
measured  by  the  Lourence  and  Pruitt  psychrometer 
system  (14),  appears  sufficient  to  give  accurate 
estimates  of  B  for  the  Bowen  ratio  technique.  The 
success  of  the  Bowen  ratio  method  for  the  estimation 
of  hourly  E  flux  rates  indicates  that  the  micro- 
meteorological  measurements  of  radiation,  temper- 
ature, and  humidity  are  representative  of  conditions 
prevailing  at  the  test  site. 

The  success  of  the  Bowen  ratio  methods  under 
poor  fetch  and  sometimes  variable  conditions 
indicates  that  fetch  requirements  are  less  strigent  for 
this  approach  as  compared  with  the  aerodynamic 
methods.  To  give  accurate  results,  the  required  height 
of  the  sensor-to-fetch-distance  ratio  appears  to  be  no 
greater  than  1  to  50  and  possibly  smaller.  On  two 
occasions,  one  over  corn  and  the  other  over  bare  soil, 
a  height-fetch  ratio  of  about  1  to  10  was  adequate. 
However,  on  October  16.  with  wind  speeds  greater 
than  500  centimeters  per  second,  a  ratio  of  1  to  75  was 
inadequate  at  one  site,  whereas  the  other  sampling 
site  with  a  ratio  of  1  to  1 50  provided  an  excellent  fit 
with  evaporation  measured  by  the  lysimeter. 
Perhaps,  during  this  last  sampling  period  a  greater 
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sampling  frequency  and  a  shorter  computational 
period  would  have  improved  the  estimate  of  E  by 
equation  3. 

On  the  basis  of  results  with  aerodynamic  estimates 
of  E  fetch  distance  appeared  to  be  from  75  to  100 
meters  for  light  and  moderate  wind  conditions,  up  to 
150  meters  or  greater  for  winds  in  excess  of  500 
centimeters  per  second  when  the  upper  profile 
measurement  height  was  1 00  centimeters  above  the 
crop  surface.  A  diabatic  correction  function  must  be 
applied  to  the  accepted  Thornthwaite  and  Holzman 
aerodynamic  expression  (29),  equation  4,  during 
nonadiabatic  conditions  in  order  to  agree  closely  with 
lysimeter  and  Bowen  ratio  E  rates.  The  Keyps  and 
Davis  empirical  correction  functions  appear  to  be 
equally  satisfactory;  however,  occasionally  a 
significant  difference  between  the  two  approaches 
occurs.  No  recommendations  can  be  made  on  the 
basis  of  the  tests  described  here  as  to  the  advantages 
of  one  correction  function  over  the  other.  Generally, 
the  aerodynamic  method  gives  better  estimates  of  E 
than  does  the  Bowen  ratio  method  only  during  early 
morning  and  late  evening  hours. 

The  error  analysis  tends  to  indicate  that  much 
uncertainty  in  the  aerodynamic  methods  results  from 
error  in  measuring  or  in  estimating  D.  This  error  is 
shown  to  increase  in  range  as  both  D  and  crop  height 
(h)  increases.  Thus,  unless  D  is  carefully  determined 
by  accurate  wind-profile  analysis,  the  aerodynamic 
methods  may  have  a  laige  error  when  applied  over 
tall  vegetation.  Thus,  aerodynamic  methods  appear 
to  require  better  conditions. 

Estimates  of  evapotranspiration  with  the  improved 
versions  of  the  combination  equation,  using  Z0 


values  based  on  wind  profile  or  crop  height,  are 
unsatisfactory  once  crop  height  exceeds  75 
centimeters.  A  ±  25  percent  uncertainty  in  Z0  results 
in  a  probable  error  of  ±  20  percent  in  E  (August  27), 
whenZ0  =  20  centimeters  as  compared  with  a  5  per- 
cent error  in  E  if  Z0  were  assumed  to  be  measured 
without  error.  Since  the  uncertainty  in  the 
determination  of  Z0  increases  with  crop  height, 
resulting  in  a  high  degree  of  uncertainty  in  potential 
evapotranspiration  estimates  over  tall  vegetation,  the 
use  of  the  log-profile  wind  function  should  be  limited 
to  bare  soil  and  to  crops  below  75  centimeters  on  flat 
terrain. 

Empirical  methods  can  be  used  to  fit  a  usuable 
wind  function  in  the  combination  equation.  If  Z 0  in 
equation  1 0  is  used  as  a  calibration  term  as  suggested 
by  Jensen  (9),  then  Z0  in  this  application  should  be 
1.0  centimeter  for  corn  as  high  as  240  centimeters. 
This  proposition  assumes  that  no  resistance  to  vapor 
transport  within  the  corn  plant  occurred  and  that  the 
potential  evapotranspiration  concept  is  valid  as 
defined.  The  empirical  coefficients  proposed  by 
Penman  for  the  wind  function  of  the  combination 
equation  appear  to  be  generally  applicable  to  a  tall 
crop  such  as  corn  in  a  humid  region.  Penman's 
version  (17)  gives  essentially  identical  results  to  Van 
Bavel's  version  (31),  using  a  Z0  of  1 .0  centimeter. 

Two  conclusions  are  possible  from  the  limited 
results  with  the  combination  equations:  ( 1 )  There  is  a 
resistance  to  water  vapor  transport  within  the  corn 
plant  despite  high  soil  moisture  conditions  or  (2)  the 
log-profile  wind  function  is  unsatisfactory  for 
application  in  the  combination  equation  where  crop 
height  exceeds  75  centimeters  on  heterogeneous 
terrain. 
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LIST  OF  SYMBOLS 


Symbol 


Explanation 


Units 


B 
D 

E 

ep 

G 
H 

Kh 

Km 

Kv 
L 

P 

R 


g 
h 
k 
In 


Bowen  ratio,  B  =  H/E 
Displacement  parameter,  D  =  d  +  zQ 
Evaporation  energy  equivalent 
Potential  evaporation 
Heat  flux  into  ground 
Heat  flux  into  air 
Eddy  diffusivity  for  heat 
Eddy  diffusivity  for  momentum 
Eddy  diffusivity  for  vapor 
Latent  heat  of  vaporization 
Atmospheric  pressure 
Net  radiation  flux 
Temperature  at  height  (z) 
Height  above  ground  surface 
Specific  heat  at  constant  pressure 
Zero  plane  displacement 
Actual  vapor  pressure  of  air 
Vapor  pressure  of  air  at  saturation 
Vapor  pressure  of  air  at  height  (z) 
Acceleration  of  gravity 
Vegetation  height 
Karman's  constant 
Natural  logarithm 
Aerodynamic  diffusion  resistance 


cm 


-2  1 

cal  cm  s 


,  -2-1 
cal  cm  s 


,      -2  -1 

cal  cm  s 


,      -2  -1 

cal  cm  s 


2  -1 

cm  s 


cm^  s  * 


cm^  s  ^ 


585  cal/g 
mbar 


i      -2  -1 

cal  cm  s 


5C,°K 


cm 


cal  gm~'  °C  * 


cm 

mbar 

mbar 

mbar 

cm  s~ 

cm 


s  cm 


-1 
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LIST  OF  SYMBOLS— Continued 


Symbol  Explanation  Units 

r$  Surface  diffusion  resistance  s  cm-' 

s  Slope  of  saturation  vapor  pressure  mbar  C~ 

curve  s  =  de/dT 

uz  Horizontal  wind  velocity  at  height  (z)  cm  s 

zQ  Roughness  length  cm 

6  Wet-bulb  depression  °C,  °K 
e  Ratio  molecular  weight  water  vapor  to  air  1 

7  Psychrometer  constant,  =  CpP/Le  mbar  C~ 
p  Density  of  moist  air                                                   g  cm" 

\Jj  Stability  profile  influence  function  1 

A  Increment 

Subscripts 

z  Property  at  height  (z)  referenced  above  ground  level. 
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APPENDIX  A— PROPAGATION  OF  ERRORS 


To  determine  the  nature  of  error  propagation  (4). 
as  error  analysis  is  sometimes  called  in  indirect 
measurement,  consider  the  quantity  Y  which  is  a 
function  of  independent  variables  Xf ,  X2, .  .  ■  Xn 
(Y  =  f[X,,  X2,  .  .  .  Xn]).  The  change  in  Y  that  is 
due  solely  to  a  change  in  X]  may  be  denoted  by 
9Y/3X,  ,  where  dY/dX\  is  the  partial  derivative  of  Y 
with  respect  to  X!  and  where  dY/dX[  represents 
the  change  in  Y  per  unit  change  in  X|  ,  while  dxj 
represents  the  actual  small  change  in  Xj .  The  total 
change  in  Y  that  is  due  to  changes  in  all  the  variables 
is  denoted  by  the  symbol  dYm  ,  where: 


dY. 


8Y 

—  dx 
3Y 


d2  dx2 


dY 

—  dxr 
3x„  r 


(1A) 

Thus,  if  dX|  .  dx2,  .  .  .  dxn  represent  the 
uncertainties  in  the  directly  measured  quantities  x  |  . 
x2  ,  .  .  .  xn ,  respectively,  and  if  all  the  uncertainties 
are  considered  additive,  then  dYm  represents  the 
maximum  uncertainty  in  Y.  This  uncertainty  in  Y 
may  be  expressed  in  terms  of  any  desired  measure  of 
dispersion,  as  long  as  the  same  measure  of  dispersion 
is  used  for  the  uncertainties  in  the  directly  measured 
quantities  also.  However,  it  is  extremely  unlikely  that 
the  actual  error  in  a  calculated  quantity  will  equal  the 
maximum  uncertainty  expressed  by  equation  1A. 

The  uncertainty  or  possible  error  in  a  calculated 
quantity  is  represented  better  by  a  probable  value 
than  by  a  maximum  value.  The  probable  uncertainty 
(dYp  )  or  error  (4)  in  a  calculated  quantity  (Y)  is 
determined  by  the  square  root  of  the  sum  of  the 
squares  of  the  individual  terms  in  relationship  to  the 
maximum  uncertainty.  The  resulting  equation  is: 


dYP  = 


dY  , 

(  )"  dx, 

9xj  1 


3Y  9 
+  (  )-  dx- 

3x2 


(—  f  dx, 
ox„ 


1/2 


(2A) 


Thus,  the  probable  uncertainty  or  error  in  Y 
depends  not  only  on  the  uncertainty  of  each  directly 


measured  quantity  but  also  on  the  partial  derivatives 
of  the  equation  that  represents  Y.  A  variable  would 
not  contribute  toward  the  probable  error  in  Y. 
regardless  of  the  value  of  dY/dX|  .  if  the  uncertainty 
(dx  1 )  for  that  variable  were  zero. The  nature  of  the 
equation  does  not  contribute  directly  to  the  error  but 
rather  modifies  the  errors  already  present  in  the 
measured  values. 

The  following  example  details  the  procedure  used 
to  evaluate  the  probable  and  the  maximum  error  in 
the  Bowen  ratio  estimate  (equation  3)  of  E. 

By  taking  the  partial  derivatives  of  equation  3 
successively  with  respect  to  Rn  ,  G,  AT,  and  Ae,  and 
by  substituting  the  resulting  expressions  in  1A  and 
2A,  one  can  obtain  the  probable  error: 


"  dR  " 


1+B 


dG-)  + 


(-RQ  +  G) 


L  (  1  +  B  )- 


B 
(  — ) 
AT  . 


-dAT- 


)  + 


(Rn+G)  B 


1  2 


(1+B)Z  Ae 


dAe2  ) 


(3A) 


and  for  the  maximum  uncertainty 


(_J_)  dR 
1+B  1 


dG  + 

KB 


-(Rn+G )  B 

 =  dAT  !  + 

(1+B)"  AT 


(Rn+G) 
(  1+B)2 


B 

Ae 


dAe 


(4A) 


where  dR    .  dG,  dAT,  and  d  represent  the  probable 

n 

error  in  measurement  of  the  variables  R  n  ,  G.A  T,  and 
A,e,  respectively. 
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On  May  2,  between  1  300  and  1 400,  the  following 
values  were  observed: 


2  . 

R     =   0.620  cal/cm~min. 

n  2 

G  =  -0.055  cal/cm   -  mm. 


AT, 


3,77  C 


10-100. 
Ae,^  ,an    =    1.507  mbar 


10-100. 
E 


0.221  cal/cm' 


The  value  of  the  probable  measurement  error  for 
each  variable  was  derived  from  table  2.  By 
substituting  this  error  in  equations  3A  and  4A,  one 
derives  the  following  error  estimates: 

=  ±  0.012  cal  cm  -nun. 
dE„  =  ±  0.018  cal  cm2  -  min. 


Thus,  E  3  1400  =  0.221    ±  0.01 2  cal/cm2  -min. 
Equations  similar  to  3A  and  4A  were  developed  from 
equations  8  and  9  by  the  procedure  outlined  above  for 
error  analysis  of  these  estimates. 


The  maximum  uncertainty  of  E  is: 

A  Au 


dE     =  ± 
m 


In 


"2-D 
-D 


2  d  Ae 


A  Ae 


Z2-X3^ 
VX3, 


dAu 


2A  Ae  Au     (Z2  -  Z 


dD 


3  (Z2  -D)  (Zj  -D) 

it 

(7A) 


Stability  Correction 


Aerodynamic  Equation 


E  =  eApk"   Ae  Au 


[  In  (Z2  -D/Z,  -D)  ]  2  (5A) 

Letting  A  =  e^p  k2  and  assuming  that  Z  is  measured 
without  error,  the  probable  error  in  estimate  of  Y  E 
is: 


dE. 


A  Au 


Z2  -D\2 
Z,  D 


dAe 


AAe 


1  n 


Z2-D 

Zj  -D 

dAu 


2 A  Ae  Au  (Z?  -Z,  ) 


n/Z2'D\  (Z,  -D)(Z2  -D) 
I  Z  ,  D  / 


dD 


J  (6A) 


4t  = 


bgAT  AZ 


Tabs(Au> 


(8A) 


The  probable  error  (di/y  for  the  estimate  of  \jj  is: 


bg  AZ 


4  ^5  Tabs  (Au)  2 
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bgAT 


^    Tabs  ,Au)^ 
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AT  AZ 


4^5(Tabs)2  (AuK 
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(9A) 
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The  maximum  uncertainty  of  is: 


d0     =  ± 
Y  m 


bg 
4  i//5 


AZ 


TabsAu' 


dT 


2pCpK    (es  "ea  >  u 


LS 


Z-D 


dZ, 


TabsAuZ 


dAZ 


The  maximum  uncertainty  of  E  p  is: 


AT  AZ 
(Tabs)2(Au)2 


dT 


abs 


-2AT  AZ 


d  Au 


(10A) 


d(Ep)m 
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dRv 


dG 
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Combination  Equation  for  Potential  ET 
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LS 
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du 


Assuming  that  S  can  be  calculated  without  error,  the 
probable  error  in  estimate  of  (E  p)  is: 
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APPENDIX    B— INSTRUMENTATION 
AND   DATA  REDUCTION 


Lysimeter 

The  weighing  lysimeter  that  was  used  to  measure  E 
directly  has  already  been  described  in  detail  (7,  8). 
The  wide  concrete  perimeter  between  the  lysimeter 
and  the  surrounding  field,  as  originally  constructed, 
has  since  been  removed,  and  crop  planting  has  been 
brought  to  the  very  edge  of  the  lysimeter.  The 
automatic  weighing  mechanism  records  weight 
changes  in  the  59,000  kilogram  mass  to  a  2.26- 
kilogram  accuracy,  which  is  equivalent  to 
approximately  0.25  millimeter  of  water  over  the 
lysimeter  surface  area.  Average  hourly  E  was 
obtained  by  averaging  the  six  consecutive  10-minute 
weights  beginning  at  60  minutes  before  each  hour  and 
ending  at  the  hour.  This  averaging  process  in  weight 
records  removed  wind-caused  irregularities. 

Wind  Profile  Measurement 
and  Data  Reduction 

The  wind  profile  measurements  at  the  lysimeter  site 
were  made  with  four  matched  rotating  cup 
anemometers  that  had  been  manufactured  by  Wong 
Laboratories.  These  were  mounted  50,  100,  150,  and 
200  centimeters  above  the  crop  surface.  The  sensors 
produced  electrical  contact  closure  for  every  0. 1  mile 
of  wind  passage,  sending  an  impulse  to  a  Sodeco 
printing  counter.  The  total  number  of  impulses  for 
each  hour  was  printed  on  the  hour,  providing  the 
average  wind  speed  in  miles  per  hour. 

Wind-profile  measurements  at  the  central 
watershed  site  were  made  during  three  of  the  seven 
runs.  During  the  July  24  run,  a  system  of  six 
Thornthwaite  rotating-cup  anemometers  was  used. 
Impulses  were  counted  on  mechanical  counters. 
These  anemometers  produce  an  impulse  for 
approximately  every  1.5  meters  of  wind  movement. 
During  the  runs  of  October  14  and  16,  a  Beckman 
and  Whitley  system  of  three  rotating-cup 
anemometers  was  used.  Impulses  were  counted  on 
mechanical  counters  that  were  photographed  on 
microfilm  every  15  minutes.  These  anemometers 
produced  an  impulse  for  approximately  every  1.5 
meters  of  wind  movement.  The  total  number  of 
counts  for  each  hour  was  recorded  and  converted  to 
average  wind  speed  in  centimeters  per  second. 

The  methods  used  to  evaluate  D  and  z0  were 
subject  to  shortcomings  and  required  sensitive 
instrumentation,     suitable    fetch,    and  neutral 
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atmospheric  conditions.  Temperature  gradients  and 
elastic  properties  of  vegetation  combined  to  make  D 
and  z  difficult  to  characterize  for  vegetated  surfaces. 
Since  temperature  gradients  are  relatively  small  over 
vegetation,  plant  geometric  properties  are  the 
greatest  influence  on  D  and  Z0 .  Therefore,  the 
parameters  can  be  expressed  as  a  function  of  crop 
height  as  shown  in  figure  9  (12). 
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Figure  9.— Roughness  length  (70)  and  zero  plane  displace- 
ment (D)  as  a  function  of  height  of  the  corn  crop  (h) 
during  the  1961  growing  season  at  Ellis  Hollow,  N.Y.  (12). 
Values  are  for  near-isothermal  conditions  and  wind  speeds  of 
100  to  400  cm/s  at  the  1-meter  level  above  the  crop  surface. 

Wind-profile  parameters  were  determined  by 
graphical  data  where  In  (Z-D)  is  visually  plotted 
against  wind  speed  and  forced  through  intercept 
Z0 values.  Since  the  data  points  were  often  too 
scattered  to  provide  a  precise  determination  of  Z0 
and  D,  these  values  were  first  estimated  on  the  basis 
of  wind-profile  parameters  measured  over  corn  of 
equal  height  (h)  at  Ithaca,  N.Y.  (figure  9).  This 
transformation  is  possible  because  of  the  similarity  in 
"Leaf  Area  Index"  (LAI)  between  the  corn  at  Ellis 
Hollow,  N.Y.,  and  the  corn  used  in  this  study 
(table  5).  Except  for  1  day  (October  14),  the 
transposed  D  and  Z0  values  provided  a  good  fit  to 
match  the  hourly  wind-profile  data.  On  October  14, 
the  D  and  Z0that  gave  the  best  straight  line  fit 
between  In  (Z-D)  and  wind  speed  were  chosen. 
Values  of  u2  andu,  were  taken  from  the  straight  line 
plot  for  each  hour  of  heights  (Z2  -  D)  and  (Z,  -  D) 
for  use  in  appropriate  equations. 


Table  5.— A  comparison  of  leaf  area  hides  (LAI)  at 
Ellis  Hollow,  N.  Y.,  and  Coshocton,  Ohio,  locations 
as  a  function  of  relative  height  (Z/h) 


z 

Z 

Z  LAI 

£  LAI 

h 

Ellis  Hollow1 

Coshocton1 

240 

1.00 

0 

0 

227 

.95 

.08 

.08 

189 

.78 

.40 

.40 

135 

.56 

1.90 

1.90 

81 

.34 

3.20 

3.20 

37 

.15 

4.00 

3.60 

0 

0 

C) 

(3) 

1  Accumulated  LAI  from  top  of  crop  downward. 

2  Base  level  of  Ellis  Hollow,  N.Y.,  corn  crop  =  4.5  LAI. 

3  Base  level  of  Coshocton,  Ohio,  corn  crop;  LAI  not  known. 


Wet-  and  Dry-Bulb  Air  Temperature 
Instrumentation 

An  aspirated  psychrometer  system  modeled  after  a 
system  developed  and  tested  at  Davis,  Calif.,  (14)  was 
used  to  measure  temperature  and  humidity  profiles. 
The  wet-bulb  air  temperature  measurement  and  the 
water  supply  were  redesigned  from  the  Davis  system. 
Changes  were  also  made  in  the  mast  and  clamping 
system,  and  minor  modifications  were  made  in  the 
sensor-holding  assembly.  Two  separate  profile 
systems  were  used:  (1)  A  10-level  profile  system  in 
midwatershed  site  and  (2)  a  5-level  profile  system  at 
the  lysimeter  site.  In  the  five  level  profile  system,  all 
measurements  were  taken  above  the  crop  canopy. 
The  psychrometer  systems  provided  measurements 
of  temperature  and  humidity  accurate  to  ±0.020°  C 
and  ±  0. 10  millibar,  respectively. 

Figure  1 0  presents  a  cutaway  diagram  of  a  dry-  and 
wet-bulb  assembly.  The  dry-bulb  thermopile  wires 
are  encased  within  a  thin-walled,  stainless-steel  tube, 
five-sixteenths  of  an  inch  in  diameter.  These  wires 
emerge  from  the  upwind  and  extend  4  centimeters 
beyond  theend  of  the  tube.  Because  the  cool  wet-bulb 
element  lies  parallel  to  this  tube  for  some  distance, 
some  cooling  effects  can  result.  Having  a  tube  of  low 
thermal  conductivity  minimizes  the  rate  at  which  heat 
can  be  conducted  along  this  tube  from  the  region 
where  the  dry-bulb  thermopiles  are  located.  A  rubber 
sleeve  shrouds  the  dry-bulb  thermopile  elements  that 
extend  from  the  stainless-steel  tube.  This  rubber 
sleeve  also  fits  snugly  over  the  end  of  the  stainless- 
steel  encasement  tube.  The  wires  are  encapsulated 
within  this  sleeve  with  silicone  rubber. 

The  dry-bulb  sensor  consists  of  a  five-junction 
thermopile  set  referenced  to  the  dry-bulb  element  of 


output  wires-l  pair  from  wet-bulb- 
depression  thermopile,  5-poir  dry^oulb 
thermopile,  wires  to  adjacent  db  s 

water  feed  tube  OD-79cm(5/l6in)- 


threaded  brass  coupling  - 
wlth "o"  ring  seal 


stainless  steel  dry-bulb 
tube  OD.-(5/l6)  in 


rubber  encapsulation 


-  dry-bulb  thermopile 
one  set  to  adjacent   dry-bulb  sensor 
one  set  to  wet  bulb  element 


Figure  10.  -Cutaway  view  of  wet-bulb  and  dry-bulb  assembly 
used  in  psychrometer  system. 

other  sensors.  An  additional  independent  five- 
junction  thermopile  set  extends  back  through  the 
mounting  plug  and  into  the  wet-bulb  element  con- 
tained on  the  same  sensor  plug  assembly.  The  dry- 
bulb  temperature  differential  system  shown  in 
figure  1 1  permitted  a  more  precise  measurement  of 
the  dry-bulb  temperature  profile  and  use  of  the 
recording  system  already  on  hand.  At  the  lower  level 
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Figure  H. — Dry-bulb  thermopile  probe  wiring  arrangement. 
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sensor,  dry-bulb  temperature  was  referenced  to  an  ice 
bath  located  at  the  base  of  the  mast,  or  the  dry-bulb 
temperature  was  referenced  to  ice  with  a  small 
thermocouple  that  had  been  encased  along  with  the 
dry-bulb  thermopile. 

The  wet-bulb  element,  which  is  shown  in  figure  1  2, 
is  not  a  conventional  cotton-wick  type,  but  a  porous 
ceramic  cup  14.6  centimeters  long. 

Water  is  supplied  to  the  porous  cup  through  a  clear 
plastic  tube  that  leads  from  a  supply  bottle  located  on 
the  main  sensor,-mounting  arm.  This  tube  is 
connected  to  a  0.79-centimeter  (five-sixteenth  of  an 


inch)  diameter  stainless-steel  tube  that  passes 
through  the  sensor  assembly  mounting  plug.  The 
porous  wet-bulb  element  is  attached  to  this  stainless- 
steel  tube  with  a  threaded-brass  connector. 

One  advantage  of  the  ceramic-cup,  wet-bulb 
element  is  a  closed  water  system.  This  system 
eliminates  part  of  the  problem  associated  with  the  wet 
cotton-wick  system,  where  a  free  water  surface  was 
required  to  supply  moisture  for  evaporation  on  the 
wick.  The  ceramic  wick  does  not  flood  when  the 
elevated  waterhead  is  subject  to  evaporation,  and 
when  the  waterhead  is  slightly  negative,  water  still 
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continues  to  move  into  the  ceramic  element  because 
the  element  acts  as  a  tensiometer  with  a  continuous 
water  column.  However,  one  must  bleed  all  the 
entrapped  air  in  the  wet-bulb  water  system  before 
starting  to  operate.  Under  most  field  conditions,  the 
water  bottle  reservoir  does  not  require  refilling  more 
often  than  every  24  hours. 

The  wet-bulb  thermopile  wires  run  inside  the 
porous  ceramic  cup  to  a  point  near  the  forward  end 
where  the  sensing  junctions  are  located.  These  wires 
leave  the  water-filled  media  through  a  hole  drilled  at 
an  angle  from  the  back  of  the  epoxy  mounting  plug 
and  through  the  side  of  the  0.79-centimeter  (five- 
sixteenth  of  an  inch)  diameter,  stainless-steel,  water- 
feed  tube.  This  exit  hole  is  later  sealed  by  refilling 
with  epoxy  after  the  sensor  wires  are  properly 
positioned.  The  five  thermopile  wires  extend  into  the 
dry-bulb  tube  that  has  been  described  earlier.  The 
two  wet-bulb  thermopile  output  lead  wires  originate 
from  the  wet-bulb  end  and  pass  from  the  sensing  plug 
with  the  thermopiles  and  follow  the  dry-bulb 
thermopile  wires  that  lead  directly  to  the  recording 
system. 

Figure  12  shows  a  scale  drawing  of  a  sensor  unit 
and  its  radiation  shielding  tubes.  A  double  shield 
arrangement  that  has  been  patterned  after  the  Davis 
design  is  used.  The  outer  shield  consists  of  an  epoxy 
fiberglass  tube  0.12  centimeter  (0.047  inch)  thick, 
while  the  inner  tube  is  a  standard  1 -inch-type  1  PVC 
pipe.  The  outer  surfaces  of  each  shield  are  covered 
with  a  layer  of  adhesive  mylar -coated  aluminum  foil. 

Air  is  drawn  past  the  dry-  and  wet-bulb  sensors  ata 
velocity  of  500  centimeters  per  second.  A  small  1  15- 
VAC  fan  capable  of  moving  1  5  cubic  feet  of  air  per 
minute  while  requiring  7  watts  of  power  does  this.  Air 
is  also  drawn  between  the  inner  and  outer  shields  ata 
rate  of  16  centimeters  per  second  through  six  small 
bleed  holes  located  on  the  back  end  of  the  inner  shield 
tube.  This  aspiration  insures  removal  of  any  heat 
loading  caused  by  radiative  or  conductive  sources. 
Figure  13  shows  complete  dry-  and  wet-bulb  unit 
along  with  all  of  its  disassembled  components. 

A  test  was  performed  to  determine  the  time- 
response  characteristics  of  the  dry-  and  wet-bulb 
elements  in  order  to  achieve  good  phase  relationship 
between  the  wet-bulb  and  the  dry-bulb  temperature 
responses  when  the  units  were  subjected  to  abrupt 
changes  in  temperature.  A  time  constant  of  20  to  25 
seconds  after  an  abrupt  change  of  7.5°  C  was 
necessary  to  achieve  62  percent  of  the  final 
temperature.  During  construction,  each  thermopile 
junction  was  checked  for  proper  insulation  and  for 
verification  that  its  response  at  a  given  temperature 
was  within  ±  10  microvolts  of  one  another.  The  final 


Figure  13. — Complete  dry-  and  wet-bulb  unit  along  with  its 
disassembled  components. 

calibration  constant  for  the  5  junction  thermopile 
assemblies  was: 

°C  =  4.876°  C/mv 

Two  profile  masts  were  outfitted  with  the 
temperature  sensors  described.  Figure  14  shows  the 
level  masts  with  the  various  sensor  spacing  used.  A 
Condorf  steel-channel  tower  10  meters  long  was  out- 
fitted with  sensor  units.  The  sensor  units  were 
mounted  with  1-inch  clamps.  Allen  wrench  set  screws 
permitted  adjustment  of  the  sensor  units  in  a  vertical 
direction.  The  dry-bulb  reference  thermopile  sets 
extended  about  2  meters  from  the  mast's  base,  which 
had  been  placed  in  an  ice  reference  bath  contained  in 
a  large,  foam-plastic  ice  pail. 

The  differential  millivolt  signals  for  temperature 
difference  and  wet-bulb  depression  were  recorded  on 
a  multichannel  millivolt  recorder  (L  &  N  Speedomax 
Type  G).  A  summary  of  these  measurements  is 
contained  in  table  6.  The  recorder  was  connected  by 
shielded-copper  leads  to  the  differential  thermopile 
system. 

Instantaneous  microvolt  readings  of  temperature 
difference  and  wet-bulb  depression  were  taken  from 
the  strip  chart  records  at  5-minute  intervals  and 
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Figure  14. — Complete  five-level,  dry-  and  wet-bulb  temperature 
mast. 


averaged  on  an  hourly  basis.  These  data  were  then 
converted  from  microvolts  to  degrees  Celeuis,  and 
average  hourly  dry-bulb  and  wet-bulb  depression 
temperatures  were  computed.  By  knowing  the  dry- 
bulb  temperature  and  wet-bulb  depression  at  each 
sampling  level,  computing  vapor  pressure  of  the  air 
for  each  hour  period  was  then  possible. 

The  dry-bulb  air  temperature  (°C)  and  vapor 
pressure  (g/m3 )  were  plotted  against  the  In  (Z-D), 
where  D  was  determined  from  wind-profile  analysis. 
A  straight  line  was  drawn  visually  to  provide  the  best 
fit  for  the  data  points  for  a  given  hour.  Values  of  e  2  , 
e  | ,  T2  ,  and  T,  were  taken  from  the  straight  line  plot 
for  each  hour  at  heights  (z2  -D)and(z.  -D)foruse 
in  the  required  computations. 


Description  of  Other  Temperature 
Measurements  and  Data  Reduction 

Other  temperature  measurements  made  at 
lysimeter  Y102C  are  summarized  in  table  6.  Briefly, 
all  absolute  temperature  data  were  recorded  on  a  16- 
channel  Leeds  and  Northrup  speedomax  G 
thermocouple  recorder  with  a  range  of  -20°  to  1 20°  F 
and  a  sensitivity  of  0.25°  F.  Dry-bulb  and  wet-bulb 


temperature  measurements  were  made  in  each  of  the 
five  aspirated  sensor  assemblies  described  in  the 
temperature  and  humidity  profile  section.  These 
measurements  provided  an  absolute  temperature 
measurement  for  the  thermopile  psychrometic 
system  and  served  as  a  check  in  comparison  with 
absolute  temperature  computed  for  each  profile  level. 

Instantaneous  temperature  data  were  taken  from 
the  strip  chart  records  at  10-minute  intervals  and 
tabulated  for  each  24-hour  period.  All  readings  were 
converted  from  degrees  Fahrenheit  to  degrees  Celsius 
and  average  hourly  values  were  computed  from  the 
10-minute  values. 

Soil  surface  temperature  was  measured  by  the 
Barnes  IT3  infrared  thermometer.  The  infrared 
thermometer  instrument  calibration  was  obtained  by 
observing  the  millivolt  output  of  the  thermometer 
sensor  while  reading  a  black-body  surface  of  known 
temperature.  The  complete  curve  for  the  working 
range  of  the  instrument  was  determined  by  repeating 
the  procedure  at  various  temperatures  in  a  water  bath 
calibration  source  that  duplicated  black-body  con- 
ditions. The  output  signal  from  the  1T3  console  was 
recorded  on  a  Royson  electronic  integrator  that  had 
an  accuracy  of  ±0.1°  C.  for  the  complete  system.  The 
millivolt  output  from  the  Barnes  infrared 
thermometer  was  averaged  for  each  hour,  and  the 
equivalent  temperature  in  °C  was  determined  from 
the  calibration  curve  of  temperature  versus  millivolt 
output. 

Description  of  Radiation  Measurement 
and  Data  Reduction 

A  summary  of  the  radiation  measurements  made  at 
lysimeter  Y102C  is  given  in  table  7.  All  the  sensors 
were  purchased  commercially,  except  for  the 
miniature  net  radiometers,  which  were  fabricated  at 
the  Water  Conservation  Laboratory  at  Phoenix, 
Ariz.  The  manufacturer's  calibration  had  to  be 
accepted  for  each  instrument  since  no  source  of 
calibration  was  available.  However,  by  comparing 
the  output  from  similar  sensors  exposed  during  a 
period  of  steady  solar  radiation,  a  field  calibration 
can  be  achieved. 

All  the  radiation  sensors  were  located  at  sites  that 
were  mostly  free  from  any  obstructions  above  the 
plane  of  the  sensing  element.  Fritschen  net 
radiometers  Nos.  432  and  567  and  the  inverted 
Eppley  No.  7583  Dl  were  all  mounted  on  the  end  of  a 
2-meter  boom  that  was  pivoted  from  a  4-inch 
diameter  mast  4  meters  high.  The  mast  was  located  on 
the  north  side  of  the  lysimeter  with  the  boom  con- 
taining the  sensors  extending  to  the  south  over  the 
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Table  6.— A  summary  of  temperature  and  wind  measurements  made  at  Coshocton,  Ohio,  summer  1969 


Observation 

Sensor 

Value 
recorded 

Recording 
system 

Sensitivity 
of  record- 
ing system 

Height 

of 
sensor 

Location 

Air  temper- 
ature (shield- 
ed and  ven- 
tilated). 

Thermopile 

Instantaneous 
mV  output  at 
5-min  inter- 
vals converted 
to  temperature 
°C  (av.  for  hour). 

L  &  N  speedo- 
max  G;  mV 
strip  chart 
recorder. 

0.005°  C 

5  levels  above 
crop  surface. 
Some  specific 
sensor  levels 
given  on  daily 
profile  data  of 
(C),  (G/m3), 
(cm/s). 

Lysimeter  Y102C. 

Wet-bulb 
temperature 
(shielded  and 
ventilated). 

-  -  do  -  - 

-  -  do  -  - 

-  -  do  -  - 

0.010°  C 

-  -  do  -  - 
(10  levels 
in  ail). 

--Do." 

Watershed  WS  101 
site  center. 

Air  temper- 
ature (shield- 
ed and  venti- 
lated). 

-  -  do  -  - 

Instantaneous 
temperature  at 
end  of  each  10- 
min  period 
(av.  for  hour). 

L  &  N 
speedomax 
G  tempera- 
ture strip 
chart  recorder. 

0.25°  F 

5  levels  above 
crop  surface. 

Lysimeter  Y102C. 

Wet-bulb 
temperature 

l  c  ri  i    I  f\ ■inn 

ventilated). 

Thermo- 
couple with 
ceramic  wick. 

-  -  do  -  - 

-  -  do  -  - 

0.25°  F 

--do-- 

--Do.  -- 

Soil 

temperature. 

Thermo- 
couple. 

--do  -- 

-  -  do  -  - 

0.25°  F 

0.05  m  below 
ground  surface. 

--Do.  -- 

Soil 

surface 

temperature. 

Barnes 
infrared 
thermometer, 
IT-3,  #493. 

Integrated 
hourly  av. 
temperature 

Royson 

electronic 

integrator. 

0.1°  C 

1.0  m  above 
soil  surface. 

--Do.-- 

Wind 

1  Wong 
cup 

anemometer. 

Av.  wind 
speed  for 
each  hour. 

Sodeco 
printing 
counter. 

1  count  for 
each  0.1- 
mi  wind. 

2.4,  2.6,  3.0, 
and  2.0  m 
above  crop 
surface. 

Watershed  WS  101 
(all  runs). 

2  Thorn- 

thwaite 

cup 

anemometers. 

-  -  do  -  - 

Mechanical 
counters 
(visual 
reading). 

1  count  for 
each  1.5-m 
wind. 

2.4,  2.6,  3.0, 
3.8,  4.6,  5.4  m 
h  =  2.4  m. 

Watershed 
WS  101 
(7/24/69). 

3  Beckman 
and  Whittley 
cup 

anemometers. 

--do  -- 

Mechanical 
counters 
(microfilm 
recording 
@  15  min. 

-  -  do  -  - 

1.8,  2.3,  3.3  m 

(10/14) 

h  =  1.3  m. 

0.5,  1.1.  1.6  m 

(10/16) 

h  =  0  cm. 

Watershed 
WS  101 
(10/14), 
(10/16). 

Wind 
direction. 

Beckman 
and  Whittley 
vane. 

Compass 
point  av. 
for  hour. 

Strip  chart 
recorder. 

±10° 

10  m  above 
soil. 

Lysimeter  Y102C. 
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Table  7.— A  summary  of  radiation  measurements  made  at  Coshocton,  Ohio,  summer  1969 


Observation 

Sensor 
and 
serial  no. 

Calibration 
constant 

Value 
recorded 

Recording 
system 

Sensitivity 
of  recording 
system 

Height 
of  sensor 

Location 

Mv/langley 

Langley/min 

M 

Net 

radiation. 

1  Fritschen 
miniature 
net  #432. 

3.12 

Integrated 
hourly  from 
continuous 
measurement. 

Royson 

electronic 

integrator. 

0.005 

1  above 

crop 

surface. 

Lysimeter 
Y102C. 

2  Fritschen 
miniature 
nef  #567. 

3.42 

-  -  do  -  - 

-  -  do  -  - 

.0005 

-  -  do  -  - 

-  -  do  -  - 

3  Fritschen 
miniature 
net  #548. 

3.43 

-  -  do  -  - 

-  -  do  -  - 

.0005 

-  -  do  -  - 

Center  of 
treated  water- 
shed WS  101. 

Incoming 
total  hemi- 
spherical 
radiation. 

Beckman  & 
Whitley 
total  ven- 
tilated. 

1.27 

-  -  do  -  - 

-  -  do  -  - 

.001 

-  -  do  -  - 

Lysimeter 
Y102C. 

Soil  heat 
flux. 

Beckman  & 
Whitley 
soil  heat 
flux  plates 
#540,  #541. 

7.7 

-  -  do  -  - 

-  -  do  -  - 

.00015 

0.05  below 
soil  surface. 

-  -  do  -  - 

Incoming 
solar 

radiation. 

Eppley  180° 
temperature 
compensated 
#7734. 

5.70 

Continuous 
trace  and 
integrated 
hourly. 

L  &  N  speedo 
max  H  and 
integrator. 

-  .0015 

4  above 

crop 

surface. 

-  -  do  -  - 

Reflected 
solar 

radiation. 

Eppley 
precision 
model  2 
#758301 

5.35 

-  -  do  -  - 

L  &  N  speedo 
max  H  and 
disc  inte- 
grator. 

-  .0015 

1  above 

crop 

surface. 

-  -  do  -  - 

(inverted). 


lysimeter.  Thus,  the  sensors  were  completely  free  of 
obstructions  except  for  the  mast,  which,  for  practical 
purposes  being  to  the  north  of  the  sensors,  had  no 
influence. 

Many  radiation  sensors  were  connected  directly  to 
electronic  integrators.  Each  integrator  system 
consisted  of  a  Royson  lectro-count  which  converted 
the  incoming  millivolt  signal  to  a  pulse  rate,  the  pulse 
rate  being  directly  proportional  to  the  magnitude  of 
the  input  signal.  This  pulse  output  from  the  electronic 
integrator  was  then  totalized  on  Sodeco  printing 
counters.  The  Sodeco  units  were  programed  to  print 
the  total  counts  for  each  hour  on  a  paper  tape  and  to 


reset  to  zero  after  printout.  Thus,  the  integrated  value 
of  radiation  from  a  sensor  could  be  obtained  by 
multiplying  the  count  output  for  each  hour  by  the 
integrator  calibration  constant  and  then  by  the  sensor 
calibration  constant. 

Two  Leeds  and  Northrup  Speedomax  H  strip  chart 
recorders  were  equipped  with  disc  integrators  and 
General  Electric  printing  demand  meters.  These  units 
had  been  assembled  by  Eppley  laboratories  for  use 
with  their  Eppley  pyranometers.  The  output  from  the 
printing  demand  meter  gave  the  total  value  of 
integrated  radiation  for  each  half-hour  period  in 
langleys. 
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APPENDIX  C-AVERAGE  HOURLY  PROFILE  DATA 


This  appendix  lists  average  hourly  profile  data  for 
each  day  on  which  observations  were  made.  Data  are 
presented  in  the  following  order: 

I.  Average  hourly  lysimeter  evaporation  and 


energy  balance  components  at  lysimeter  and 
midwatershed  sites. 

2.  Smoothed  profile  data  from  log-linear  plots  of 
temperature,  humidity,  and  wind  at  lysimeter  and 
midwatershed  sites. 


31 


CO 

uu 
I— 
< 

cc 

x 
3 
-j 

LL 


<T 

LL 
X 


O 
Z 
< 


7* 


< 

Of. 

(_ 

Q. 
< 
> 
U 


3 
IS. 


> 

— 

c 
I 


< 

cc 
Lt 
> 
<3 


C] 

3 


21 

cl 


O 
X; 

uj 


u 
I  i 
CM 
o 
—t 
>l 

a 

d 

H 

ILt 

3 

3 


•a 


I 

i 

uu 

< 

c 


t 

\  CM 

I  z: 

!  U 

iS  — « 

I  _i  r» 

!  <  —  o  o 

I  O  I-  o  o 

I      <!0  O 

!  I  o  •  • 

iO  o 

;  co 


o  o 
o  o 


o  o  o 


! 

I 

O'O 
oo 


o'o  o!o  o  o  o!o  O  O  0| 
o|o  OiO  O  O  O  O  CIO  o 
OIO  OiO  O  O  O'O  o;o  o! 


o  o 


O  O^O  o 


!  a 

■  Cl 
X 

a 
o 


sCiO  O  O  — i  CM 
«— !  O  Ol-^-  OM^ 

colo  o!o  o 
o1  •  • 

a;jo  o 


o  o 


in  i  CO 

— 41  CM 


o  o  o oo  o  o  o  o o 


. — i  [     i— ii  f- 

rv||  f\)  ogi  — ' 


ro  o  Mo  oio 

O!"4  -0  I  fVl  O'O 
— ii— '  OiO  O'O 


O  OjO  OiO  OIO  0:0  OIO  C 


p  o  I  co  cm 
o  — » j  o  o  m  —• 

L0|0  O  -4  *t 
O  O 


i/iiO  n 


— l|vO 


— i '  in  in  |  rvj 
in  |  cm  ^o 
M  cm  — i|  o 


o;— t  in  oc  inio 
t  |  m  «o  Um  o  o 

co  sC  roj— i  o'o 


o  o  — •  -< 


— I   —I  O 


O'O  C;0  O 


I  o 

!  <3 

I  a: 


2:  <fiir. 


— '  o  — • 

J-  o  r- 


— -!  — 1  — I 


I 


|  2. 

1  |  a:!o 
cm  it 

U 

j  "v.  j<j 
|  _J  —  M 

<  roio 
lUH  • 

!  CNO 

I  •  i 

i 

|  LU 

;  o  io 

|  Z*  —i  co 

;   <    C\||  — I 


OIO  o 

I  I 


rv1:— 1 

O.'lO 


rvj  a>  Oii^ 
■0  o  sC  ;  in 


oo  j  o  .m 

^  (V1  c 


o  o  o  olooio 


o  cc  I  U  CC  I 


o  O  OiCJ  O  O  Oi 
!        I  I    I  I 


O  o|  rv  Njvn  in;  m  in:  -T  <+•  ro  — i*  — j  o!  O  Oj 

O  OiO  O  O  OjOO  O  0;0  O'O  o'o  oi 


o 


ol  o 
I  I  I 


cr  i  oo 

oi  o 

Oi  —i 


O  O  Oi  O  O  O  O  O  OiO  o  o  o 

ill    III    til  Ml 


rn'cT^  vC'  o  cm;  •— <  so  co  co  r*-  r-iCT1  r~ 
m;(\J  O!  vj-  sT;  (T  O  ^  <V'  0s  -J-1  «t 
— i  i  r^>  <rl>j-ro|ro  cm  i  •— 1  oi  O  — 4 !  — '  rH 


<  n\  o 

CO 


> 

(  0 


LL  — '!  O 
2T  LUi  • 

LU  _j!o 


o  c! 


o!  o 


o  o 
o  o 
o  o 


o  o 

'1 1 

I 

o!  rvi 
in 

O  -4 


O!  O 


0;0 

I  I 


0  o 

1  I 


0  oio  o1  o  o  o 

1  Ml     I  I 


f»-  pn  r\j  cm  j  o  O  |  o  o :  <_> 
^0  c\)!  — *  —I'm  o i  O  Oi O 
rvj|  — -  — i'o  Oio  o|o 


0  o 

1  I 


0  o 

1  I 


o 


1  LU  Q- 1  IT 

LU ' 
!  LL  Q-| 

S  J 

I  i-  ol 


s0 


0  c 

1  I 


O  O'O  oio  o  j  o 

III  i 


r-  co  0s  c 


rv  m  <ri  in  o  s  cc 


o  o^— 4 
— >  cv  cv 


o 


CM  I 


LU 

u 
< 

— 
3 
CO 


c 

OCT 


> 

<l 

3 


Xi 
3 


<i 

LUi 

X 


00| 

II I 

— ij  ! 

LL  O: 
0O         — i 


!  +1 

2t;  o! 

j  id 


Ui 


1  ; 

II: 

<l: 

x| 

coi 
i — < 

H 

>i 

•* 

_J 

rv  j 

v! 

Ll'; 

3 

1— 

_a 

i— i 

Oi 

Ujl 

LLi 

i — j 

U4 

.X: 

<r' 

1— 

»— ii 

1—' 

uj 

ir\ 

<r! 

U-l 

>! 

uJ 

<r 

_i 

_J 

X 

«j 

u_ 

LU 

LU 

i— 

of 

c' 

_J 

LLI 

OS 

CO 

II 

LLi 

1 

II 

oo 

LU 

aj 

— f 

LU 

zj 

DC. 

—i 

r  i 

•J- 

1 

2rj 

3j 

cc\ 

"I 

M 

i.u| 

a 

Oi 

-5. 

<Ti 

ri 

i 

Zli 

-—>; 

3 

ol 

a| 

3 

oo. 

3 

J 

0: 

<ri 

0 

col 

U| 

_J: 

Oi 

d 

1 — 

1— ll 

d 

sj 

I 

c| 

•—1 

•! 

11 

Of; 

<J1 

oo 

— J1 

a 

32 


< 

cc 


U4 


< 

LL 

J 

I— ) 

c 
<y 

c 


q 


a  3 


2 

C 


< 

Cr 
C 
Q 

< 


►4 

4 


c 
I 

C 


c 


i 

in 


11 

< 


2 

1 1 

i  <m 


UJ 

Z 

Q 
CL 
V 

O 

o 


q  O 
H  o 

<!  O 
°1 


!  i 


0  o 

ot  o 

01  o 


o]  o 
Oi  o 
o|  o 


c  o  o  o 
o  o  o  o 

O  O!  O  O 


o  o|oo 
o  o'  o  o 


o  o  o  o  o  o  o  o  o  o 


O  O;  O  O 

o  otoo 


o  o  o 
o  oj  o 
o  o  o 


oqoci 


-Oj  o 
«-«(  o 
ooi  o 
o  • 
o 


^C!  m 


O  O  Oi  o 


LO  mi  m  lo|  O 


<f  m!  O  m 


>  o  o  o  o 


o 

Q 


r- 

o 

a 

f\j 

o 

o 

Oi 

o 

o 

o 

o 

o 

° 

-  lO 

o  Jo 

<—  ini  o 

<r  ooj  o 

o 


0|  LO 

(V  o 

ol  o 


mi  r- 
r-i  r- 

r\i! 


x 


mi  m  oi     r*-!  r—  oi  fv  o  o  m  o 


2 


_i  _j       m  ,£>:  an  f-v  CM       O  Ol  O 


o  o 

i 


Oj  O  O  O  O1  O  Oi  O  O  O  Ol  O  Oj  o 


o 

oi 

°i 
j 

o 


3 

.5) 
X 


3 

5 


a 

Mi 

a 


a 


N  N  (\j 

i  | 

,  i 

:  I 


z! 

Q 
»-l 

H 

«1 


a 


a 
-I 

a 

j 


a 

3 


>—i 

O: 

o 

rr\\ 

ooj 

r— 1  CO 

<r 

o 

Loi 

d 

z 

cc 

»±i 

*i 

m 

cu 

\0  f\l 

o 

o 

-1 

-I 

o 

oj 

Oi 

i — i 

mi 

<r 

m|  »v 

o 

o 

-1 

— < 

H 

2. 

• 

• 

• 

m 

1 

Oil  o 

oi 

c 

oj 

O 

o 

o 

o! 

o 

oi 

o 

oi  o 

o 

o 

oi 

o 

o 

2" 

M 

1 

ij 

1 

1 

I 

i 

1 . 

1 

1 1 

! 

1 

i 

a 

•i 

+: 

llj 

u 

i 

1 

zi 

a) 

r*-i 

LT 

rsji 

m 

<r 

>— lj 

r~- 

r- 

>o 

<r 

— * 

m 

+ 

LO 

— 1 

—J 

===== 

H 

o 

o 

o 

o 

o 

— Jj 

m'  rv 

o 

•—i 

LU 

a 

a 

< 

ro 

o 

o 

o 

O 

o 

° 

o 

Oi 

o 

oj 

o 

ol  o 

o 

o 

o 

o 

o 

LO 

_l 

CO 

o 

• 

O 

o 

a 

o 

Ol  o 

o 

Oi 

o 

O 

o 

C>!  O 

o 

C- >; 

o 

o 

1 

in 

_j 

1 

1 ! 

1 

T 

1 

I 

I 

i ; 

1 

1 1  1 

1 

T 

1 

II 

*-j 

<t 

Li.1 

i 

xi 

i 

u 

rxj 

— < 

i 

Is- 1  m 

-rv 

sT 

vD 

r~ 

f\J 

lt.i 

■C 

f\)i 

CO 

n> 

* 

3 

— i 

m 

LOj 

—1 

m 

sC 

Oi 

<r 

(\i 

r— 1 

mi  so 

o 

m 

rvj 

m 

o 

cy: 

X 

O 

< 

f\l 

— i 

, — 1 

oj 

#— i 

— < 

o 

o 

,  1 

— <|  C3 

— 1 

- 

LU 

3 

a 

_J 

H 

-1 

a 

<r 

X 

o 

oj 

o 

oi 

o 

o 

o 

oi 

o 

o 

o 

O;  O 

o 

O 

oi 

o 

o 

lU 

LU 

iJ 

UJ 

CO 

1 

7' 

r 

1  |  1 

<j 

H 

»— 1 

d 

U 

> 

oj 

to 

< 

UJ 

<r 

o 

1*- 

LTV 

o 

o 

Cf) 

LO 

LfM  O 

o 

rv 

00i 

CO 

CJ 

>- 

LU 

-1 

Cf' 

— < 

o 

o 

o 

o 

vol 

CO 

<r 

Oi 

o 

<r 

<— 4 

r— 1 

o 

_J 

LLl 

LU 

o 

o 

O 

o 

o 

O'  o 

cvj 

m 

i— (i  — ' 

. — 1 

O 

oi 

o 

o, 

1 

UJ 

• 

• 

5; 

LU 

LU 

o 

a' 

O 

O 

o 

o 

o 

a' 

o 

C 

o 

o;  o 

o 

o 

o 

O 

o 

a 

—1 

1 

i 

1 

i ! 

1 

1 

!  ! 

1 

1 

1 

1 !  i 

1 

1 

1  i 

I 

a; 

cO 

II; 

LU 

►—1 

c 

i 

i 

lO 

m! 

LO 

LJ 

LU 

3" 

i 

—1 

LU 

g 

LU 

LO, 

r» 

cd 

<^ 

o 

rsii 

m 

LO 

cc 

C7- 

o 

r-4 

rv 

cO 

04 

I 


33 


3 


M 


I 

PJ 


<  w 

U  h- 

< 


oo  oo 

(Nil 

in  vO 


h-j m  oojcr  on<n  cr>  rP  r-  r-  —uro  trr^  pj> 
<f-jpj  -o:>f  r*-!o  O'fM  oo;rp  o;oo  m!ro 
c?v  — <  m,sD  r*- 1 r*-  <o! *o     ;  pj  o  1     -j-  <f  ^  j 


O  OIO  Ojr-i  — H  i- i  — i;— i  — i!-*  OfO  0,0  O: 


— 

O  GO 
Q  O 

a  i 
o  I 

n  !° 
CJ  — »•  o 

i—  LOl  o 

V—  — •  • 

<r  oo! o 
-  o! 

a  : 

<  ! 


rp  O 

o!  lp 

O!  O 


on!  — *  O;  LP 
CO'  ■— i  LP  i  no 


aojO  o  co  o>s  <M;^  miO  oj 

0*i  ;  CP  r»-  irp  r> !  LP  OO  Oi 

. — i !  — *     4 ;  •— i  i— » ! r-i  OiO  O  O  O i 


oio  oio  O'O  o  o  0:0  o  o  o:o  oio  o, 


ooj  LP  rpl  o 

rpl  00  r— ti  f\j 
LP  CO;  —* 


rvj 

r— I 

fM 

O  OI  O  O  v— 


o  lp  o  co  r-  ■  o  o:r^  co  o 
co  •  lp  ro|o  r- ,  lp  cvi!o  rp  i  o 
—i|  o  cr  :r-      pj  0|0 


M  H  i— l|r-l  O'O  0,0 


o; 
ol 
o 

• 

o! 


I*- 

ro 

CO 

f\l 

vC 

rp 

o 

—1 

LP 

rP 

LP 

o 

— i 

LPI  CU 

o 

ro 

LP 

CO 

, — 1 

<_J 

LP 

CO' 

■*■ 

-i 

<^ 

c 

pj 

LP 

X 

^- 

LP 

PJ 

o 

. — 1 

1 

ed  o  o 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

(_> 

o 

o 

o 

o 

o 

r\j 

1 

1 

1 

1 

1 

>: 

. — i 

O 

<r 

r» 

<r 

LP 

sT 

PJ 

u^ 

r» 

— " 

PJ 

_j 

o 

o 

rv 

ro 

-p 

PJ 

PJ 

PI 

rp 

ro 

pi 

o 

Oi  CJ 

o 

< 

mi  c 

Ol  o 

o  o 

O 

O 

o 

o 

O 

o 

O 

O 

Q 

o;  o 

o 

u 

i 

u 

o 

o 

1 ' 

O'  o 
1  !  1 

o 
1 

O 
1 

o 
1 

o 
1 

O 

1 

o 
1 

o 
I 

o 
t 

o 
1 

o 
1 

o 
i 

o 

o 

UJ 

O 

vo 

oj  o 

— i 

i— < 

PJ 

<r 

PJ 

IT 

LP,  vC 

LP 

jC 

CO 

rP 

z 

lP 

CO. 

Pj;  sO 

m 

O 

ro 

a. 

o 

r- 

o 

ppi 

LP'  -O 

r- 

< 

o 

o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O: 

< 

X 

O 

o  o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

Ol 

CC 

1 

1 

1 

1 

> 

an 

i 

o 

Ol  PvJ 

00 

00,  IT, 

r~ 

PJ 

r- 

PI 

PJ 

CC  PJ 

■Jj 

oi 

0£ 

— i 

o 

r- 

p- 

ir\ 

PJ 

ro 

rn  "p 

o 

r- 

1 — 1 

1 — 1 

LP 

ol 

UJ 

o 

pj 

PJ 

LP 

r- 

CC 

ar 

CO 

00 

O  LP 

PJ! 

O 

oj 

z 

UJ 

UJ 

—J 

o 

o 

1 

1 

o 
1 

O 

1 

o 
I 

o 

1 

o 
1 

o 
I 

o 
1 

o 
1 

o|  O 

1  i  1 

1 

c> 
1 

o 
1 

Oj 

a 

i  a  o 
z »-! 

I  LL  CZ\  LP 
I  LU 
i  X.  CL 
S  ! 
'  Luj 


>Gl  r-  CC  0^  Oi  — '  P* 

I  <|  r-H  r— J 


rp  >a-|Lr,  vOjp-  cc      o,-h  pj; 

—t  r— <  i  t— i      \  — i  ■— i  <— i  pj  pj  rv; 


LU 


— ' 


— i 

I 

A 

<3 
•j 

q 
H 

H 

a 


3 


crj 

rj 

d 
5; 

o 
< 

UJ 

< 


H 
U 

u 


< 

CC. 

u 

od 
< 

LU 

I 

&! 

H  oJ 


a 

S3 
or 


< 

a 


I 

CM 

2l  , 
U 

v.  — 
_j  r- 
<  — 

a  h 


O  O,  W  D!  m  iri  vO  vO  — i  o>d-  O  o 
Ma^-HiOOON^^  -^--IP   O  i  rO  !N  CO  r- 

co  av  — *  lo;  r-  oci  ao  cx>.      -c'  o  co;  r—  n0:  lo  lo1 


O  O  O  O'  — I  — i  :  — i  — ■>       O  O,  O  O  O  O;  o  o 


LU  sO'  o 
I  2T— *|0 
|  O  c/J|  O 
lo.  o  . 

^  a:  o 
!  O  I 

!° 

z  o 

lo  -Jo 
! «  m  o 

<  ml  o 

i—  oi 

l«  ! 


! 

— j  ro 

Z  <f|  vO 

-Ho 

*  2  * 

I  ad  o 

tN)  j  | 
5. 

*s     I  <r 

-j  o 

<   CV  O 

u  —I  • 

Oi  o 


O .  ro  !_n  to  _ij  P«- 
O  O  rO!  r—  vC.  -C 

o|  C  O!  O  Hi  H 


f-j  lp  ro!  lo  r<^ :  lo  oj;  lo.  lo  o 
n\  ro  ccj  r>-  mi  r—  m|  p°  — o 
— i  —  -<|  O  Oj  O  O;  O  O,  O 


C !  O  0:  O  O]  O  Oi  o  ol  o  oi  o  o  o  o  o 


0|  ro  CO,  O  roi  LO 

c  s  m  OD  M  O1 


o  o 


1*—  LO,  ro.  rO  fM  rO  sO  LO  LP 

— *  lo.  — o  io  o  lo  m  ro  o 
m  3oio  n     h  >j-  m;^  n  -  o  o 


ol 
o; 
ol 


O'  O  O  O  oj  — <  — i  — i  O  O IO  O;  c  oi  O  Ol 


>f  N  -f !  <I  -f  M 

i.o  — j  vtv  -r  r\ji  cp- 

O'  O  oi  Oj  io:  -O 


O  o  o  o  o 


vC      ^c>T|      u>  a  ^  n  cr  M 

LO;  LO  vO  LO-  0:^7)  O  — i  f-  ^-  <LT 
CO  rX3  CO  Oj        r\:        O  O  O  O' 


O 


O  O  O  O  O  o  o  o  o  o  O 

III 


I 


-ij  rti  r\i  >u  o,  h  r-i  a;  to  r—  co:  — •  i—  r\j  ,_i  — i 

Oi  o  o;  O  ■— I       ro  rv  O'  — i  o!  — i  U."1  O  O  O  O 

o  o  oi  O  Oj  o  o  o  o  o  o:  o  o  o  o  o  oj 

o  o  o:  o  oj  o  o  o  o  o  oi  o  O  O  O.  O  Oj 

I     !  I  1  I     I     I     I  !  Ill  I     I  I 


ill    III     III    I  I 


I   I  I 


O  Q 

^  n 

LL-  CC  \  lO 
LUj 

lu  oJ 

3C  i 


-C  r-  x 


CT  Oi        OJ  rr>  ^j-   LO  >£,  P-  CO  CT  0;  — *  fvj 


LU 

a 
< 

LU 

OlI 

Z3 
'^1 

51 

a 

ILi 

> 

<; 
< 

X 

o 
_l 

LU 

^1 

LU 

< 


1—1 

t/i 


UJ 

1 

U 

I  a* 

LO 

so! 

r— 

ro 

o 

CO 

vO, 

*J  LO 

O; 

r— ( 

o 

~*  oci  <f  vO  O  O  — i  lO:  rrv 

m 

LO 

o 

r~- 

o 

<: 

f\j  O 

o 

O! 

O  ojl  o  oj  o  O  rv  oi  c_> 

O 

1 — 1 

oj 

ai 

— j 

o  oi  o  o  o  o  o  oi  O 

i— 

<r 

o 

Q 

Oj 

o 

o 

O; 

o 

LU 

OJ 

cc 

xj  o 

T 

1 1 

i  1 1  i  i ;    r    i ;  i 

> 

-J  o 

P«-i 

o 

o 

<r  c!  m  cc          c  o>  co 

~Oi 

cc 

(\|J 

O 

ol 

>i 

Ql 

H  ° 

LOi 

o 

LO! 

OJ  OJj        LO:  rr-.  r*^-  sj-  O;  vO 

LO. 

— i 

o 

o 

_J 

UJ 

—  o 

o 

o 

oi 

OJ  OJJ  sO  Loi  CD  CO   sj-  oi  ~* 

o 

l—l 

o 

i! 

35 

JJ 

— 1  ° 

o 

o 

Oj 

o  o!  o  o  o  o  o  ol  o 

o 

c 

Oi  o 

Oj 

a 

LU 


LU 


lil 


cc 

>< 
3 


< 

LU 
X 


o; 

cj 
i 

Z' 

3 
1 

oo 

Uj 


U| 

O 
> 

d 

LO 

I 

ad 

LU 

H 

LU 

5; 
■— i 


sO 


a 


z  | 

It  i 

CM  | 

I  a:  i 
K  J 

-J  M 

i  <  —:  o> 
jo  hj^ 

i  <, 

II  O!  • 

j  I O 

i  00  I 

1 1— 
Z  ~l 

:  UJ  \0  O 

Iz  w  o 

iO  wo 

j  0.  31 


v0        Hi  n0 

o  <i-  r-'ro 
m  h-  o  co 


cc'-f  coir-  o  ao  cm  cm  o  o  lo 

mi>o  O; <J-  ro  (M  o  co  cm!nO  m 

LOisO  r»-;<o  loco  olr-  if»  >t  -4-ivi- 

•    •    •    •    •    •         •    •  • 


O  O  M 

I 

!  i 

CO  s0  LO  s0 

ml  o»  4*  i*- 

O  O  — i 


O  OIO  O 


o  o 


c 
u 


o  o 


Z  lO 

lo-Jo 
I  —  lo!  o 


o  o 


>3- 

o  m  o 


<  ooj  o 
a 

!  «a  ! 


r-  O 


Oj  23 
CM  I  —I 


OicO 


o  ro 

C  1  LO 

— "  o 


ro|  O  O  O  O 

o|o  OIO  o 
o;o  oio  ol 


o  olo  o  o  o  o  o;o  oio  o  o  o 


LOiLO 

co  o 

0:  r-l 


ro-  cm 
h! 


cc:  r- 

O !  CO 


co  to 
r-ir- 

LO  I  CM 


|  i 
0:0  O'O  O: 

oi c  o:o  o 
o! o  oio  o! 


o  o  o 


. — i  f . — »  ^iO  O'O  o;o  o'c 


z  >* 

«— i  — «j  o 

5-  Z|  • 

I  0£0 

fM  !  I 
X 

<->  1 


O  LO|  CO 
m  LO  CD 
— *  ro  LO 


o:r- 

o!o 
r—i  ao 


Ci  rr, 
cc;  r- 


>t:  ro 

P-1 


r\j I  co 
o*»j  vc 
ro  j  i—i 


i       i  > 

CJM  CM  NO  CO: 

CM1  <3  P*"  !  LO  %0| 

O  O  OiO  Oi 


o  o  o|o  oio  olo  o;o  oio 

I 


I  <  ro;  o 
IO-  • 
o  o 

I  I  ! 

I 

!  ^  | 

:  o  cr 

1  Z  — J  LO. 
'  <  CM|  O 


■c  co  — ( 
-•  o  o 

Ol  o  o 


Oi  o  o 


cMj  r\j  p— 
o  o  o 


LTj  r— 
■— li  CM 

olo 


mj  cr 
roj  f—i 

Ol  o 


>r  o 

o:  O 


I 

o:o 
o!  o 


o|o  OIO  O' 
Ml    ll    I  ! 


ro  r~  a;|  4j  cu 

O1  O  o:  O  oj 
o!  o  o; O  O; 


i  <  x:  o 

I  m 
I  > 

;  o  — I  o 

CX  — li  o 
LU  — I  O 

!  z  uj;  • 

i  UJ  _|i  o 

I  I 

i  ol 
I  O  qJ 

!  z  — 

lu  a 

LU 
U_  Q. 


0  o 

1  i  I 

i 

crj  co 
p-ir- 


o;  o 
I  i  I 


rot  >1" 

LO1  O 
r\l  r-4 


O1  o 

I  I 


v0  ro 
ro 


0  o 

1  I 


LO  nJ- 

o  o 


ol  o  o 
I  I  I 


Oi  CM  OIO 
Oj  i-H  CO  I  CC 
I  —*  CM  I  CM 


or 


0  o 

1  :  I 


CM  |  f\, 

o  o 

LO,  LO 


o1  o 

i !  i 

I 

col  O 
lo  r— 

LO  v0 


O  O 

1 1  I 
a 

f— •  X 

s0;  CM 


O'O 
I  j 


O  CO 
CO'  >0 
CM'  r-* 


Ol  O  O  o  o: 


CO i  sO  >D|  «0  C 
-£>  LO  LOi  LO  O' 

H  o  o  o  o 


I 


o  o  o  o 

I  I  1 1  I 


0  o 

1  I 


OlO 

i !  i 


0  o 

1  I 


olo 

I  :  I 


o  o  oi  o  o 
IjllM 


cc  a*lo  f-i.(\j  ro.  <f  lo:  vo  r-!co 


CT!  O  -h  rv  ro 

fij  CM  CM!  CM  CM 


LU 

< 

U- 

31 

LO 


CC, 
LL. 

>-i 
< 

Xi 
Oi 

-I 

LL. 


O1  O  O  O  O.  -H 


sTi  ro  O  Oi 
O'  CM!  O  s0: 
>— i  _ <  -h  oi 


< 

3  o; 
J  < 
a.  -£! 

h-j  ac; 
< 
_jl 

oi 


X 
-J 

C 

1 


o1 


<; 
o 


o 


Oi 

< 

LU 
>i 
<| 

51 


<: 

| 


Si 

o| 

Oi 


oo 

o 

LU 

ro| 

Zj 

< 

1— (i 

O 

<i 

U 

I 

C! 

1 

oi 

r- 

d 

zi 

21 

Z| 

H 

H 

» 

51 

ill 

o; 

23 

Oj 

O 

+! 

-! 

i— < 

LU 

Oj 

^_ 

-J 

'^0 

—1 

" <j 

1 ; 

LO 

_l 

1 

in 

< 

1 

CO 

*1 

i 

a 

*| 

J 

•i 

x; 

zi 

XJ 

E3 

H 

J 

a 

UJ 

LL, 

H 

LU 

Sj 

<i 

\— 

»-* 

o 

<i 

aj 

>- 

LU 

<s 

—I 

_J 

LL 

s: 

LU 

i— i 

LU 

or 

C5 

_J 

LU 

CD, 

z 

li 

1 

00 

III 

LO 

lU 

Z1 

O 

-I 

uj 

z 

ry 

a: 

— i 

— 1 

CM 

36 


< 

DC 
X 

LL 


< 

LL 
I 

J 
h-i 

d 
d 


4 

3 

1 


g 
S 

LLi 


51 

H 
-4 


LL 

C 

< 

a 

LL) 
< 


c 


c 

X 

o' 

L_ 
Lj 


I 

c 


tz- 
ar 

LL 
h- 

LL 


> 


4r 


O 


li 

M 
< 

a 


I 

oj 


<  — ! 
o  (~, 

<! 
I  Oj 


I  ° 

!  o 


i-i  in 

K  — 

<  00 

•—  C? 

a 
a: 


CO 

N0 

o- 

CM 

4- 

**  in 

■C 

• 

• 

o 

o 

o 

O 

o 

a 

in 

m 

c 

o 

m 

o 

3 

o 

o 

• 

• 

o 

C 

o 

o 

a 

o 

m 

CO 

o 

o 

cc 

o 

o 

o 

_l 

• 

• 

• 

o 

3 

o 

o 

m  roi  CO  Oj  r- 
O  CO'  4"  CJ»!  4^  O 

r*-  sO  o-  ooi  co 


m  co!  O  co 
<\i  4-j r-  <■* 

in  ;  m  m 


Oj  O  Oi  O  O!  o  O!  o  ol  o  o 


4" 

m 


o  in 


h-  cc|  — i 


I 


*  m  Ni  o 


m|  -4-  oi  o 
-o  4-;  o  o  c 


O  0|  O  OI  O  Oi  o  o!  o  o;  o 
c  o  o!  o  o:  o  o:  o  oi  o  oi  o  oi  o 


co  r*-i  «o 

LO  4"j  O 

oj  oj;  m 


mi  c\j  >0;  m      m  O:  o 
Ol  — •  col  lo  oj  4"  Oi  O 
4-  cvj;  m  oj!  o  o,  o 


o  o  o  o;  o  oj  o  oi  o 


r- 

M 

O 

4" 

LO 

?"•» 

I— lj 

lO 

4 

m 

mi  cm 

cc  c 

z 

4" 

m 

mi 

m 

>c 

m 

rMj 

m 

o 

in 

ir\j  n 

° 

o 

-< 

m; 

CM 

1 — 1 

r— 1 

o 

i — i 

— <:  r— 1 

z 

• 

• 

• 

• 

•j 

• 

•J 

• 

• 

• 

•  * 

1 

o 

o 

o 

o 

o 

o 

o 

O: 

o 

o 

o 

Ol  o 

f\J 

! 

lj 

1 

1 

1  1 

5. 

O 

j 

\ 

lj 

oj 

u 

o 

o 

o 

Oi 

u 

oj 

o 

CJj  CJ 

Oi  <U> 

-pi 

o 

O 

o 

o 

o 

o 

Ol 

o 

O 

o 

o 

o 

o  o 

< 

m 

o 

a 

o 

o 

o 

o 

oi 

o 

o 

o 

O'  o  o  o 

o 

• 

«i 

• 

•1 

• 

• 

• 

• 

• 

• 

•  • 

1 

5 

Oi 

o 

o 

o 

o 

o 

Oj 

o 

°" 

o 

O;  O 

LU 

o 


< 

CD 

> 

O 
CC 
LU 


o 

CM 

o 

a* 

-Oi 

m 

c- 

m 

CM 

m 

m 

CM 

r- 

4- 

o 

i— i 

Q 

o 

o 

r-4 

o 

o 

XiOOOOOOiOOOO.  oo 
Ml     l|l     III  II 


H  >0 

H  LO 

— ■  o 

m!  • 

1° 

I 

a 


oo  <r  in 

O  •— 'i  i— ' 

oj  ^vJ;  OJ 

o  o;  o 


N  M  ^  O  U  Cl  IA  N  U  >C  *i  -U  O,  >C 

•— •  •— i  CO  LO  >J1  OJ  CCj  — H  0^;  lO  LO:  LT\  LO  lO 

-hi— ioo-^-^— i  •— •  a  o  o!  o  oi  o 
«  •  H 


oj  o  Ol 
I   I  l! 


qo  o|o  do 
hi    HI    I  I 


o;  o  c  o 
l>l    I  \ 


r*~  cc  c>j  o  H  cm 
i  ■— i  i— i  ■— < 


^0  h"  co 

I— I  r- 1!  r— 1  r— f. 


o 
< 

LL 

z> 

s: 
a 
a: 

LL 

> 

< 


<     !  O 


i  X 

:  — J 

| 

UJ 

I- 

< 


a    j  — 


o 
oo 


LU 


OO 

> 


5" 
O 
a: 


2 

a 


i  < 


i  < 

i  o 

i  < 

I  cc 

LU 

!  > 


X 

1  i— 

s  o 

!  3T 

<i 

LL 

;  a: 

i  _ 1 

— J 

i 

1  rv 

\  cc 

< 

<l 

LU 

j  _J 

:  o 

X 

!  a 

1  z 

_ J 

]  oo 

»— 1 

i  ^ 

r-> 

'  _J 

1  o 

oo 

< 

u 

; « 

1! 

!  u 

I  — 1 

:  — i 

O 

i  ° 

I  'i 

— . 

o 

I— 

m 

< 

1—1 

— 

o 

w 

z: 

— 

«b 

+ 

!  1! 

C£ 

<I 

+ 

oo 

_J 

LU 

:  o 

a 

— 1 

00 

,  

1 

LO 

II 

< 

X 

j  CO 

i  o 

» 

-J 

X 

o 

a 

c 

LL 

r- 

LU 

< 

'  (— 

o 

< 

UJ 

LU 

< 

!  _i 

X 

I  cc 

LL 

LU 

LU 

i  H 

CC 

— i 

LU 

cn 

:  m 

Hp 

00 

'  II 

i  oo 

!  O 

LU 

a: 

O0 

:  OC 

i  vO 

OJ 

4" 

37 


OO, 
UJ 


21 
I 

(X! 


i  ! 

1  ^ 
1  -j 

Pi 

i 

I*1 

— i>t 

0- 

<*■ 

X 

O 

X! 

° 

r-;ro 

O  j-< 

m 

tvj 

3  i 

<  in 

in  ih- 

CO 

o,<> 

lf\ 

— 1 

i  1 

o:  • 

* !  * 

• 

• 

* 1  * 

• 

• 

• 

• 

LU 

lo 

o;o 

O 

0 

OO 

0 

O 

0 

0 

O 

1 

i  00 

■  t- 

j 

< 

z 

LU 

Lu 

m|o 

ro  in 

conn 

m 

r\l 

r\l 

in 

0! 

X 

! 

1  z 

w  ,— 1 

cm \r- 

rH 

m 

m 

(\l 

*—> 

O: 

0 

oo;o 

00 

OiO 

O 

0 

O 

0 

O: 

a. 

!  2: 

cc  jo 

o;o 

O 

O 

010 

O 

0 

O 

O 

Oj 

OO 


9 


a 
< 
> 

OJ 

1 

OO 

d 

X 
uJ 

a: 

LL 
> 


Ci 
X 

o 


o 

Q! 

Ij 
00 


u 
I 

Oj 
>-| 
Lb' 


LU 

M 

LU' 

*— *; 
OO 
>-! 


i— I 

VJ 
r-4 
I 

I 

UJ 

H 
< 

Qj 


a 


d 

i  !z 

X 

m:c\! 

p— < 

r» 

cc 

in 

0 

in 

zj  ! 

1  ;a 

<!- 

cm 

in 

CO 

sO 

CM 

0 

OlCVJ 

0- 

m 

CM 

0 

*  i  * 

• 

• 

• 

• 

• 

■ 

< 

00 

0 

oio 

0 

0 

O 

0 

O 

0 

o 


I  1-4  Qj 

o  1 
!<  1 

(X 


rsj  jr— 

c 

<n 

r» 

C\i 

C1 

0 

!  2" 

1 

m  r- 

0 

LP 

m 

. — 1 

r-i  m 

0 

1 

— jo 

(— 

r«o 

— 1 

. — ; 

0 

O'  0 

0 

Is. 

!  1 

cx.,0  OtO  O 

O 

0 

0 

0 

0  0 

0 

j  <NJ 

1  !  1 

1 

2 

i 

1 

0  <_> 

<_> 

l_j 

CJ 

U!<J 

LJ 

i  -J 

— .'O 

O'O 

0 

O 

c 

0 

0 

o;  O 

O 

'  < 

rnio 

OIO 

0 

C 

0 

0 

0 

0  0: 0 

O 

!  1 

0,0 

OIO 

0 

LJ 

0 

0 

0 

0  u 

O 

I 

j  LU 
i  O 

ir» 

inj— ' 

cr 

LO 

. — 1 

m 

r— 1 

0" 

o;r^ 

cr 

1— 1 

1— 1 

rOi  O 

IS 

CMjO 

1— 1 

0 

0 

O 

1— < 

— 1  f— * 

00 

0 

O 

0 

0 

O 

0 

O  O 

0 

i 

! 1 

1 

1 

l 

|  > 

i 

1  a 

«— .!  r-l 

— < 

O 

m 

cc 

(Ml  vD 

! 

vO'tn 

0 

CT 

-T 

m 

0 

p-4  in 

in- 

:  LU 

—  <— 1 

oj  0 

0 

I  Z 

\^ 

o!o 

0 

O 

O 

0 

0 

0 

00 

0 

I    I    I  I 


O' 

0 

D| 

z 

i 

LU 

a  k 

u«j 

LU 

0-! 

»— 1 

Ltl 

(- 

Oj 

I  II   111  I  I 


38 


Q 

-£ 

O 


> 

o 

1— 

o 

r- 

o 

rC| 

o 

•J- 

o 

o 

LO 

m 

PC 

(NJ 

p^ 

LO 

i— < 

LU 

t— ( 

00 

O 

1— 1 

PC. 

or, 

rn 

0s 

co 

r— 

r- 

lT 

C3 

■— I 

f—i 

(NJ 

OJ 

m 

rn 

rc 

PO 

-n 

PC 

(NJ 

— 1 

i— ■ 

1 

2: 

O 

' — 

— ■ 

LU 

o 

21 

0£ 

—* 

— 1 

CO 

CC 

r- 

h- 

(N) 

CO 

LO 

>o 

v0 

—i 

O 

o 

LO 

(NJ 

PC 

o 

LO 

>0 

r- 

X 

(NJ 

LO 

\0 

Q* 

LO 

LO 

r>" 

LO 

sO 

pc 

CO 

i — i 

1— 

™j 

, — i 

f\i 

(\i 

CM 

(NJ 

M 

!NJ 

>— ( 

■— i 

• 

<I 

o 

Cl 

LU 

Q_ 

2! 

II 

LU 

r~ 

c 

o 

LO 

CG 

>T 

LO 

o 

sT  O 

1  1 

rn 

—^ 

\ 

o 

rH 

C\J 

CO 

vD 

O 

CO 

to 

>o 

LO 

0 

^-^ 

PC 

PO 

X 

LL 

5- 

r- 4 

o 

L_> 

O 

LO 

LO 

to 

LO 

yC 

O 

~o 

r*- 

f~ 

CO 

CO 

CO 

c& 

X 

l—l 

O 

• 

LU 

CO 

*™ 

2*. 

X 

h- 

X 

>■ 

CJ 

o 

H- 

0_ 

-J 

o 

o 

Cu 

Q 

o 

ct 

c 

r  > 

(  ) 

(Nj 

1 

PV 

t\ 

o 

m 

fc 

PC 

cc 

2 

• — i 

C\J 

CT 

>£) 

CO 

ro 

O 

»r 

(NJ 

>T 

pc 

o 

(NJ 

o 

LU 

o 

X 

LO 

LO 

LO 

LO 

•O 

p*- 

r*- 

r~- 

CO 

CO 

CT1- 

o> 

cr 

CO 

CO 

CO 

X 

< 

t 

— ^ 

o 

lo 

i 
1 

cc 

L_ ) 

O 

O 

-J 

o 

o 

■C 

fC| 

PC 

r— 

LO 

p— ( 

•J- 

PC, 

CO 

LO 

(NJ 

o 

II 

• 

o 

o 

cc 

^O 

PC 

3^ 

>3" 

rn 

LP 

o 

>t 

CO 

p*- 

o 

5. 

L/) 

(_J 

o 

■  ■ 

LT 

co 

E — ) 

f\l 

fC1 

<r 

LO 

<T 

(NJ 

o 

CO 

H— 

3£ 

1 

C\J 

CVJ 

OJ 

(NJ 

(NJ 

(NJ 

CNI 

(NJ 

Z" 

II 

LL 

<NJ 

"— 

LJ 

o 

5" 

Q 

<J 

1— t 

LU 

LU 

k- 

> 

CC 

rx 

<T 

oj 

z> 

C 

LU 

o 

LO 

o 

3D 

m 

m 

O 

OJ 

r- 

CO 

LTi 

LO 

(Nl 

LO 

o 

(NJ 

OJ 

L/l 

(— 

< 

. — 1 

O 

LO 

r~ 

pc 

o 

cc 

CO 

LP. 

o 

r- 

CO 

r- 

LO 

<r 

■d 

LU 

i— * 

or 

LU 

—J 

00 

in 

LU 

CO 

X 

& 

r\l 

r- 

OJ 

<r 

LO 

r» 

C0 

CC 

cc 

—i 

CO 

LO 

O 

Cl 

. — i 

CM 

oj 

r\) 

PM 

f\i 

(NJ 

(NJ 

OJ 

(NJ 

(N 

r— i 

— i 

—4 

LL 

C£ 

LL 

O 

o 

LU 

i 

LU 

a 

—4 

a" 

t— 

I— 

• 

a 

UJ 

LU 

h- 

O 

s: 

X 

o 

>— * 

< 

a 

O 

II 

LU 

CO 

o 

o 

I 

z: 

LU 

O 

y- 

_J 

LU 

DC 

LO 

CO 

o 

CO 

LO 

sC 

CO 

o 

i—l 

(N) 

X 

rxl 

O 

LL 

rv 

(NJ 

(NJ 

u 

LU 

Q. 

5l 

to 

M 

LL 

>- 

□ 

39 


□ 
< 


** 

>- 

i 

f— 

o  o 

LLi  r— J 

Q 

uO 

SN 

■5: 

J~ 

*" 

LJ 

UJ 

0 

CC 

t— *  A 

— ' 

L- 

«-r 
<J. 

CY 

Li_< 

CI 

1  1  ] 

pO 

"5  ( — 1 

^  w 

IX 

r*  i 

i — i 

>- 

C  3 

I  I 

i 

o 

- 

L — 

I — i 

m  CD 

?~ 

;  o 

.._ 

LiJ 

o 

-w- 

X 

' y") 

• 

1 

CJ 

| 

[5 

CD 

<^ 

O 

/ — \ 

I  j 

1 

L  3 

u_ 

.  (Nj 

C 

<T 

|JJ 

CH 

< 

CM 

C 

LU 

(XI 

*~  0 

h- 

< 

LL 

i— * 

\ 

cc 

_J 

oo 

to 

UJ 

a 

UL 

ex. 

:> 

c 

UJ 

1 

LiJ 

cc 

h- 

t— 

a 

Uj 

UJ 

? 

o 

< 

O 

u 

a 

Q  j — 

X 

> 

Z  1— 

y- 

 1 

c 

UL 

CJ 

uj  a. 

5 

2L 

1—  O 

LOr<^ojcT'Ln'J>r\jO(NJLr\'>nos--^'-^Lno<J-Nt 
_4  _i  _i  — 1  — 4(\i(NiCN|rM^(\ir\j(Mf\j^jrvi 


r^iv-r^-coovosovC7sc>c^oNCN00craoxih- 


in-iOs+oo^Kinf\io^^Hf\j^CTN 
r^f^ocCT-oocT^cro— 1  m  (\  - '  h  o  a1  x  r 


oomoocoiri^"nirlcr'tcra>oOxOX\r 
X'xi(3,''-|(M'NJ'-|-lNfri^in>j-vriAj^(jo 


in  J3  s  x  a  o  - i(NjrOsj-ir,  >o  h-  00  CT 


O  — *  CM 

f\J   !N1  f\J 


o 

o 


< 
II 

UJ 

5: 

UJ 
DC 

D 


X 

UJ 
X 


o 


X 


a 

0 

o 


2.' 


o 

o 


11 

G 


2: 
-J 

O 
r- 1 
• 

o 
II 

o 


40 


> 

o 

O 

o 

O 

O 

r- 

r— 

LO 

CM 

LT, 

00 

CC 

O 

O 

r- 

LO 

ro 

LLi 

CM 

o 

O 

■— i 

—* 

03 

LO 

cc 

rO 

(<fl 

r*» 

o 

i— i 

CM 

O 

00 

oO 

— 1 

CM 

fM 

CM 

CM 

rO 

ro 

>f 

ro 

ro 

CM 

CM 

CM 

—4 

« 

">v 

5" 

O 

X 

~" 

#» 

LU 

O 

C£ 

in 

ro 

IT, 

LO 

<T 

O 

LO 

CO 

O 

O 

o 

o 

o 

LO 

cr 

CM 

ro 

u 

X 

ro 

ro 

O 

r- 

r- 

rO 

<T 

CM 

vC 

>c 

i— i 

i — i 

LO 

r- 

co 

LO 

k- 

— i 

—4 

— 

i — i 

— 1 

*- • 

CM 

pr 

o"i 

cO 

ro, 

CM 

•— 1 

i — i 

—4 

• 

< 

LO 

r- 

lu 

C 

5l 

ii 

LU 

PO 

o 

LO 

Is- 

r- 

o 

LP 

IK- 

O 

CO 

o 

-T 

O 

LO 

CM 

O 

CM 

h— 

\ 

LO 

ao 

m 

00 

CM 

CC 

O 

O 

o 

CM 

CT 

i— 1 

r1^ 

CM 

X 

Li- 

5: 

CM 

o 

c 

r— 

r- 

r- 

CO 

CT 

O 

c 

O 

o 

c- 

cr 

cr 

cr 

cr 

cr 

cr 

O 

* 

CM 

(M 

CM 

CM 

— < 

i— i 

t— 1 

i-H 

i—i 

< 

—> 

— * 

• 

LU 

OO 

X 

h- 

X 

>- 

C 

X 

i— 

-J 

C_J 

CJ 

CU 

o 

LP- 

O 

c_J 

cr 

o 

CM 

LO 

LO 

O 

o 

o 

cr 

o 

— •* 

CM 

o 

cr 

LO 

p-- 

i— ( 

X 

LO 

00 

c 

••C 

00 

-o 

cr 

—4 

ro 

ro 

r<^ 

< 

o 

_J 

o 

LU 

o 

X 

r- 

cr 

1  1 

•— ^ 

— 4 

~ 1 

o 

cr 

cr 

cr 

o 

cr- 

cr 

cr 

cr 

X 

CM 

CM 

CM 

CM 

CM 

i — i 

. — i 

CM 

— 1 

OO 

• 

oo 

1 

cu 

UJ 

<r 

(—J 

D 

o 

_J 

o 

lo 

LP 

IT 

LO 

o- 

C 

OJ 

CM 

cr 

— ^ 

r— 

o 

LO 

o 

ro 

vf 

I! 

• 

LTv 

LO 

— - 

>g" 

LO 

— t 

cr 

>3" 

CC 

lO 

sD 

Q 

LO 

CM 

CM 

cr 

ro 

>o 

o 

CM 

00 

uO 

u 

~~~~ 

■ 

1 

— i 

co 

U*i 

o 

cc 

cr 

cr 

CJ 

— =J 

O 

CO 

! 

1 

uj 

CM 

CM 

CM 

C\J 

CM 

t\l 

CM 

<M 

CI 

C*l 

ro 

CM 

CM 

IM 

^M 

II 

LU 

f\J 

cr 

LU 

o 

-o 

21 

c 

<r 

— i 

LU 

LU 

> 

\ 

C£ 

CC 

< 

X 

Q 

LU 

CM 

t— 

o 

r- 

ro 

r- 

CM 

(M 

vO 

C 

LO 

ro 

r— 

CM 

•0 

o 

O 

LO 

LO 

00 

h— 

<r 

LO 

ro 

o 

r— 

ro 

o 

>i- 

CJ" 

LO 

LO 

CM 

cr 

—J 

ry-\ 

< 

LU 

LU 

51 

—I 

CO 

>u 

Ui 

r— ! 

ro 

LO 

co 

cr 

cr 

O 

O 

f— t 

00 

r- 

r- 

~o 

s: 

O 

CL 

Oj 

CM 

(\i 

CM 

CM 

CM 

CM 

CM 

CM 

co 

CO 

ro 

ro 

CO 

CM 

CM 

CM 

CM 

LL. 

QL 

5" 

o 

O 

LU 

I 

LU 

o 

o 

rx 

h- 

• 

CL 

LL' 

LU 

h- 

CM 

51 

k- 

X 

o 

1—4 

< 

Q 

o 

II 

LU 

a 

o 

> 

■z 

UJ 

o 

1— 

LU 

IT 

r- 

CO 

O 

CM 

CO 

LO 

r- 

00 

cr 

O 

—1 

CM 

X 

CJ 

LU 

CM 

CM 

CM 

(_J 

LU 

a 

2t 

oo 

LU 

o 

»— < 

X 


a 


•> 

> 

.... 

O 

o 

Is- 

X 

s0 

X 

-J- 

X\ 

X 

r~ 

~0 

1—1 

LU 

rn 

<f 

LO 

X 

X 

X 

rn 

O 

X 

a 

!/1 

—1 

—< 

—i 

o 

I 

»* 

c 

LU 

z 

o 

QC 

X 

(M 

1 

O 

— 1 

•— < 

O 

CM 

X 

sO 

LTi 

X 

o 

r- 

*H 

C_) 

CM 

X 

in 

<f 

m 

X* 

X 

o 

X 

—J 

X 

O 

co 

r— 

co 

r-t 

• 

< 

o 

a: 

j- 

UJ 

CM 

a 

2l 

— . 

II 

LU 

m 

h- 

O 

O 

o 

CM 

o 

s0 

in 

CM 

x 

<x 

in 

O 

o 

<M 

o 

o 

o 

h- 

*>. 

in 

a; 

i—i 

in 

X 

CM 

X' 

rn 

X> 

_i 

a 

Is- 

r» 

X 

i— i 

X 

CM 

o 

X 

LL 

21 

m 

O 

C 

O 

sO 

^0 

o 

«o 

r- 

r- 

X 

r- 

O 

Is- 

r- 

Is- 

r*- 

i— i 

C 

• 

LU 

uo 

2. 

X 

t— 

X 

> 

'O 

O 

o 

l— 

a 

_J 

o 

C 

a 

o 

X 

cc 

O 

!»— 1 

o 

CM 

r*- 

oc 

o 

X 

in 

v0 

IT 

X 

X 

X 

fM 

r- 

m 

O 

(X 

5 

in 

o 

LO 

<r 

a 

CM 

r- 

>t 

o 

<i 

u 

X> 

CM 

CJ 

Lw 

o 

X 

>o 

<J 

x 

r~ 

«o 

Is- 

r~ 

sC 

Is- 

X' 

sO 

r- 

Is- 

r- 

z 

X 

< 

_ 

CO 

• 

_i 

u 

o 

l 

u 

Cj 

u 

o 

< 

<M 

o 

o 

_J 

O 

X 

rn 

in 

XI 

ir 

CM 

x 

m 

X 

o 

n 

X 

in 

Q 

CM 

a 

II 

• 

Is- 

rxj 

in 

_« 

r- 

r- 

■c 

X 

1 — 1 

in 

CM 

o 

CM 

m 

o 

o 

2; 

m 

m 

m 

(_ 

l_> 

X 

cc 

;r 

cr 

o 

CM 

<r 

m 

-£) 

•T 

>r 

rn 

CM 

CNJ 

>— 

a. 

1 

u 

t~i 

•—i 

— i 

CM 

CM 

CM 

CM 

CM 

CM 

i\ 

CNJ 

z 

!l 

LL 

(M 

III 

o 

-O 

s: 

O 

< 

dJ 

LU 

1— 

>- 

a 

CC 

< 

D 

10 

Q 

LU 

cm 

h- 

o 

Is- 

ac 

o 

f— 4 

CM 

r- 

in 

CM 

O 

o 

cc 

X 

un 

o 

CO 

h- 

< 

X 

m 

X 

x 

m 

in 

o 

<—< 

CM 

in 

<M 

<r 

i-J 

<r 

o 

< 

LU 

\ 

CC. 

CM 

LU 

5" 

_J 

l/> 

r- 

LU 

ao 

cr- 

CJ 

ro 

o 

r- 

sO 

•<r 

<3" 

rn 

CM 

i 

a 

i — i 

. — i 

i-H 

CM 

CM 

(M 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

<M 

CM 

CM 

CM 

x 

a 

s 

LL 

o 

C 

LL 

1 

uu 

o 

o 

or 

t— 

h- 

• 

Ql 

LU 

LU 

t- 

n 

X 

—j 

c 

o. 

O 

o 

n 

LU 

a 

a 

c 

X 

> 

z 

LU 

o 

t~ 

Lb 

in 

O 

r» 

00 

O 

O 

r— 1 

CM 

m 

O 

Is- 

X 

& 

O 

—> 

r\i 

X 

c 

LU 

CM 

CM 

CM 

u 

LU 

CL 

a: 

on 

LL 

o 

42 


Q 

z: 
»— i 

a 
< 


•> 

> 

O 

u; 

LO 

LO 

f\l 

0 

co 

LO 

LO 

CM 

r- 

•4- 

>— 1 

CM 

<J- 

<— 1 

UJ 

m 

CM 

>t 

LO 

0 

—1 

CO 

rn 

•— ( 

>— < 

cr, 

O 

m 

LO 

-O 

*o 

CO 

O 

00 

.—1 

— -4 

r— 1 

—1 

.—1 

—1 

2. 

ZD 

O 

X 

•» 

O 

UJ 

0 

or 

t— t 

[f\ 

O 

LO 

ro 

O 

CM 

CM 

cr 

CO 

ro 

cj 

r-< 

CM 

<T 

O 

c 

r»- 

co 

X 

X 

10 

r- 

4- 

CM 

m 

<r 

CM 

ro 

I- 

1 — 1 

• 

< 

O 

cr. 

■4" 

UJ 

CM 

Q. 

— 

II 

UJ 

ro 

K 

51 

O 

LP 

N 

O 

CM 

T— < 

LO 

r- 

X 

<t 

O 

LO 

>o 

O 

ro 

t— 

^. 

LO 

O 

0 

0 

"Tj 

•4- 

cr 

>t 

O 

X 

r» 

CM 

X 

U_ 

5- 

ro 

(J 

O 

0> 

cr 

cr 

cr 

O 

0 

O 

O 

0 

O 

O 

O 

O 

cr 

cr 

cr 

cr 

O 

>— 1 

O 

■ 

LU 

00 

X 

1— 

X 

> 

O 

0 

r-i 

t— 

Ci 

_l 

« t 

c_> 

a 

Q. 

O 

in 

c_ 

0 

<M 

co 

r- 

cm 

LO 

LP 

LO 

—J 

CM 

t— 1 

ro 

ro 

m 

CM 

r- 

ro 

O 

h- 

ST 

LO. 

O 

— 1 

r- 

C 

m 

0 

. — 1 

cr 

X 

X 

X 

<. 

0 

_> 

CM 

O 

UJ 

0 

X 

cr 

cr 

0 

O 

0 

r  -1 

— 1 

O 

0 

O 

c 

cr 

cr 

cr 

cr 

0 

z 

X 

^—1 

< 

00 

• 

_ I 

O 

00 

2. 

1 

O 

a  \ 

c2 

< 

CM 

0 

O 

_j 

O 

LO 

O 

0 

0 

O 

O 

<T 

-O 

r- 

O 

LO 

0 

*t 

0 

II 

• 

to 

>o 

00 

-O 

<1 

cr 

X 

ro 

ro 

cr 

-3" 

<f 

cr 

r- 

r- 

Cr 

O 

5L 

00 

LO 

c; 

O 

— 

CM 

CM 

m 

<! 

X 

0 

, — 1 

m 

sT 

«T 

<r 

CM 

■—1 

f—< 

cr 

X 

h- 

a1 

1 

#— 4 

— i 

« 

—1 

1— 1 

CM 

(M 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

— i 

II 

u_ 

CM 

,  . 

LU 

O 

5T 

O 

< 

1— 1 

UJ 

UJ 

1— 

> 

\ 

a. 

CI 

< 

r» 

ZD 

lu 

r\J 

y— 

O 

ro 

0 

LO 

O 

a 

LA 

LP, 

X 

r- 

m 

LO 

<! 

O  f\J 

O 

lO 

O 

0^ 

t- 

< 

lo 

CT 

r- 

cc 

X 

>t 

«t 

0 

-0 

r- 

CM 

>r 

CM 

r- 

CT 

< 

UJ 

s 

cm 

LU 

—J 

00 

uu 

UL 

.— * 

!M 

en 

X 

O 

CM 

<*) 

-T 

<r 

<r 

IO 

CM 

cr 

X 

CJ 

t— 1 

a 

CM 

CM 

CM 

(XI 

CM 

Cxj 

CM 

CM 

CM 

CM 

<— 1 

1— 1 

UL 

cr: 

LU 

O 

c 

LU' 

1 

UJ 

O 

c 

a 

t— 

t— 

• 

a. 

UJ 

UJ 

h- 

0 

X 

X 

CM 

c 

< 

c 

II 

LU 

00 

Q 

a 

0 

X 

>■ 

LU 

O 

K 

-J 

LU 

cx 

O 

oc 

0 

CM 

m 

vT 

LO 

O 

r- 

X 

cr 

O 

1— 1 

CM 

X 

0 

LU 

CM 

CM 

CM 

CM 

u 

UJ 

Q. 

2_ 

i. 

to 

M 

LU 

1— 

O 

43 


o 

' — i 

s. 

a 
z 
< 


> 

O 

O 

LP 

ro 

CO 

CT- 

00 

co 

p\ 

—1 

PJ 

LP 

LP 

r- 

<=— < 

LL1 

fVI 

r\j 

7— 4 

0 

CT 

0 

IT' 

CO 

0 

Pj 

r — ) 

sC 

vO 

<) 

(/I 

OJ 

fVJ 

M 

C\J 

CNJ 

rv 

C\l 

CN 

C\l 

rH 

. — t 

. — i 

— - 

27 

x 

' — ' 

— 

O 

UJ 

X 

0 

"Si. 

» — 1 

in 

r— 

LP 

0 

ro 

O 

a 

(NJ 

LP 

LP 

LP 

0 

0 

~ 

IS 

~r 

r- 

h- 

(Ni 

<f 

sT 

LP 

lP 

r- 

cc 

C\i 

(NJ 

•p 

1— 

f— 1 

1 — 1 

—1 

f\l 

00 

rv 

• 

< 

0 

or 

ro 

UJ 

—< 

*r  _ 

"™" 

IU 

V- 

2" 

0 

-0 

0 

r*- 

-4" 

00 

lP 

O 

r— 

i— 

LP 

co 

<T 

(Nl 

<T 

. — 1 

O 

CT 

X 

LL 

2L 

CNI 

X 

O 

0 

LP 

LP 

LP 

LP 

LP 

LP 

LP 

P 

LP 

LP 

LP 

rP 

• 

LU 

L/5 

2. 

X 

^~ 

O 

-J 

O 

\  1 

Qu 

2. 

1 — 1 

LP 

CC 

M 

O 

iP 

r\i 

CO 

(NJ 

X 

LP 

CJ 

LP 

rv 

O 

ac 

1 — 

2: 

LP 

. — 1 

LP 

c~ 

■J" 

CO 

CL 

■p 

rn 

1 — 1 

1 

O 

_J 

Q 

UJ 

C 

X 

ip 

lO 

If 

IP 

LP 

•p> 

LP 

LP 

LP 

'P. 

lP 

z 

X 

y2 

—J 

LJ 

w  1 

1 

lj 

X 

u  \ 

LJ 

(J 

r— 1 

CT 

C 

_j 

LP 

O 

0 

r*» 

CC 

CM 

PJ 

00 

ro 

C 

X 

X 

II 

• 

in 

CO 

0 

LO 

1 

LP 

LP 

LP 

1 

CC 

nO 

(NJ 

O 

L/" ' 

LP 

LJ 

' 

LP 

CO 

CC 

X 

CO 

cr^ 

O 

O 

CO 

ro 

• 

X 

1 

LJ 

F— J 

1  ' 

II 

U_ 

(Nl 

U_' 

O 

2: 

<T 

LLj 

U- 

> 

•s 

rx 

..X 

< 

— 1 

X) 

0 

LU 

t— 

0 

or 

c 

LP 

LO 

<r 

CO 

O 

LP 

O 

O 

O 

O 

1— 

<r 

LP 

c 

<r 

r» 

o> 

ro 

h- 

CO 

PJ 

X- 

LP 

X 

<r 

UJ 

> 

cc 

=— 1 

LU 

— 1 

LO 

0 

JU 

X 

X 

cc 

CO' 

O 

0 

LT 

X 

r>- 

<T) 

0 

~4 

□l 

» — 1 

•—1 

— H 

LL 

a: 

5: 

LL 

0 

a 

LU 

1 

LU 

O 

0 

a: 

t— 

1— 

• 

(X 

LU 

LU 

I— 

X 

2" 

I— 

X 

a 

< 

O 

11 

LU 

'/I 

LJ 

Q 

c 

t— < 

X 

> 

2T 

UJ 

0 

1— 

Lb 

cx 

sC 

r- 

CO 

0s 

0 

lP 

X 

O 

X 

c 

LU 

(NJ 

u 

LU 

a 

2: 

2_ 

on 

u_ 

O 

44 


> 

I— 

o 

X 


—  a 

o  o 

LU  .-I 

OO 


LJ 


DC 

O 

0 

CM 

UTS 

AO 

lt\ 

O 

>d- 

<— < 

m 

in 

ZD 

-5 

cr- 

0 

& 

c 

<f 

— 1 

0 

cr- 

1— 

C\l 

r\l 

m 

CM 

f\J 

<\1 

— 1 

— 1 

• 

< 

LU 
Q. 

UJ 


LL 

o 


00 


1— 

X 

w 

C 

_l 

a 

z: 

c; 

^- 

< 

0 

LU 

2! 

X 

00 

_J 

! 

O 

O 

—1 

2: 

CD 

O 

1 

LL 

C\l 

CT< 

0 

>£> 

< 

1— 

> 

<i 

O 

a 

LU 

~* 

1— 

aJ 

"X 

_J 

OO 

G 

LL 

a. 

O 

LU 

1 

at 

1— 

Q_ 

LU 

LU 

2" 

1— 

O 

< 

LU 

OO 

X 

> 

1— 

_) 

O 

00 

27  O 
\  O 

O 


_3 
X 


oo 

CO 

< 


o 
o 


LU 

cc 

O 
<  —1 

Cl 
LU 
CL 


a 

Q  D 
Z  — 

LU  £X 
LU 

u_  a 
s. 

—  u_ 

o 


r-  r*-  lo  r\j  'm->j-r\jr<*itncc 


^(NJcch-^ror-vCLom'Xicr 


2: 
o 


o 
o 


a 
< 

t— 1 
11 

00 


LU 
cr 
ZD 
00 

LU 

5. 

LL 

C 


X 


LU 
X 


X 

LD 

LU 
X 

a 
0 

a: 
o 


o 

o 
o 

II 

Q 


5_ 


O 

o 


o 


45 


Q 

Z 

i 

Jt 

! 

z 

< 

; 

>■ 

i 

*■•  o 

1 

O  LO 

i 

UU  CO 

I/O 

*"* 

X 

3-  ! 

ID 

O 

X 

"** 

i 

CJ  ! 

UJ 

i 

z  o 

cc 

— 1  in 

^~ 

i 

<I 

CC 

1  I  1 

: 

UJ 

i  

r~ 

i 

ST  O 

u_ 

i.  |P»! 

1 — 1 

O  : 

CD 

U~ 

1 — 

X 

i 

' 

(  J 

I — 1 

N«*  O 

QC 

h— 

<l 

o 

ZD  rsj 

uu 

CD 

x 

4C 

•X. 

I 

r  f\ 
W  ) 

1 

L-l 

1 

I  1 

<~ 1 

/ — 1 

«J 

O 

'UJ 

JL. 

1 — 1 

LJ~ 

Lr) 

lj> 

LL 

j 

o 

xj 

<c 

UJ  j 

CC 

< 

 I 

^" 

X>  ! 

o 

f\! 

IO 

^  uo 

LU 

Q£  rsj 

■J 

k/i 

[ 

UJ 

i_ 

0  ! 

LL 

im 

O 

1 

1 1  1 

rV 
□C 

ix 

CL 

P 

UJ 

I 

O 

LU 

< 

UJ 

a 

I 

!< 

t- 

O 

a 

uj  (a. 

2: 

LO 

nL 

I 

cocoaooooO'^0N00>-^—<oacaooooooOf-i 
^vO<OvOif>vOoo(\if^Mr^-toca5ooooirLn 

■-»(   »«*  (\l  t-"4  1— I 


ooir\H>j-OHt\j(\j(N(Mif\fn(\irno(N(\j 


^OOC<\)^OOLOun<\lcC'^"'rvOr\JfNJ<\l,^3LO 

^0f<^^oc^c^olf^p,-  >otocr'-^o''-'>J_ur\(\jr^- 
hpJh-*hm«njm  r\jc\ir\]r\jr\jr\jc\jc\irvjr\i 


l^fniv-OOvfs-'^^-'-l-|,-|^^'-'^vC(\JOO 

aocococr'Ov<\j>a-Ln  o^inmrnrnw^HO 
•^^^^r^rM<\jrsir^(NrvjrM<Njr\jr\if\jr\j(\j 


o 


X 

o 


in»CNocp,OHf\jm<i-fn^ir-flcc,0'-'N 


• 

LU 

X 

O 

a  ] 

O 

a 

uo 

a: 

m 

u 

0 

< 

1 

0 

2.  1 

m 

t\j 

0 

II 

• 

O  ! 

lo 

uo 

1— 

Z 

II 

uj 
2: 

c 

UJ 

cc 

_j 

00 

i 

< 

LU 

2: 

27 

0 

LL 

0 

O 

0 

•  i 

h- 

LO 

X 

— -t  : 

O 

II  [ 

LU 

°,j 

X 

M  1 

46 


3S 
a 
< 


> 

o 

o 

ITi 

n 

rO 

CM 

CNl 

rn 

in 

n 

CNJ 

4 

r- 1 

— 4 

fM 

4" 

UJ 

ro 

4" 

a 

»— i 

X 

—j 

— i 

X 

O 

in 

■C. 

m 

4" 

Q 

L/0 

— * 

1 

—i 

i— i 

■H 

r: 

X 

■k 

LU 

z 

O 

5". 

CC 

m 

o 

c 

IT 

O 

4 

m 

CM 

CM 

<C 

•4- 

>C 

m 

PO 

CM 

— J 

U 

=> 

<M 

o 

O 

cc 

o* 

DO 

00 

IT 

h- 

(M 

rr 

<* 

CM 

ro 

— i 

• 

< 

cc 

UJ 
CL 

UJ 


O 

oo 
I— 
O 
_j 

Qu 

CC 
< 


I 

o 


2. 

u 

CC 
LL 


< 

Q 

UJ 


LL 

jo 

iCC 
Q. 

|o 

iUJ 
X 

;u 
lac 
></, 


ro 


X 

fx 

r>- 

r\j 

m 

3 

CM 

CM 

in 

end 

up 

CM 

— ( 

X 

>^ 

i — ^ 

m 

*o 

ro 

rn 

O 

TO 

cc 

X 

CT 

cr 

cr 

y 

0" 

cr 

0^ 

cr 

sO 

49 

r- 

a 

M 

X 

> 

c 

1— 

M 

Q 

c 

o 

CM 

CM 

m 

X 

r- 

o 

o 

cr 

vT 

(M 

CM 

ir 

o 

LP 

4" 

1— 

2T 

in 

CM 

4- 

a 

cr 

r- 

LT 

r- 

>£) 

X 

ro 

ro 

o 

-J 

N 

o 

X 

CO 

cc 

X 

cr 

cr 

cr 

o 

cr 

cr 

O 

cr 

sO 

r^ 

X 

—I 

•—I 

• 

a 

00 

u 

CD 

< 

o 

in 

m 

in 

>o 

N 

o 

o 

X 

r- 

O 

m 

in 

CM 

o 

(M 

in 

— < 

CNI 

m 

*o 

r- 

CM 

c- 

CJ> 

X 

o 

ro 

cr 

o 

X 

ro 

o 

CM 

in 

X 

O 

r— 1 

—i 

CM 

ro 

m 

— « 

cu 

r» 

>0 

LO 

O 

CM 

CM 

CNJ 

CNJ 

CM 

CM 

cr 

o 

LU 

CC 

r-  ■ 

CL 

ro 

i— 

O 

r- 

CC' 

4- 

X 

m 

N 

O 

p» 

ro 

c 

CM 

rv 

r- 

4 

h- 

rn 

< 

in 

1— 1 

4" 

o 

CM 

r- 

m 

CM 

4" 

vO 

CM 

< 

n. 

CC 

<M 

Q0 

LU 

o 

o 

in 

00 

cr 

o 

CM 

ro 

ro 

ro 

— 

X 

m 

a. 

.—I 

— 

CM 

C\) 

CM 

r\ 

CM 

CM 

CM 

o 

UJ 

1  1 

LL1 

H 

m 

UJ 

cc 

UJ 

< 

a 

1— ■ 

a 

C 

o 

< 

z 

3 

LU 

cc 

r— 

00 

cr 

c 

i— i 

<M 

ro 

in 

s0 

r» 

X 

O 

o 

4" 

CM 


X 


• 

UJ 

X 

u 

Q. 

o 

O 

n 

cc 

ro 

u 

o 

z. 

< 

2l 

in 

U 

CM 

o 

II 

• 

o 

m 

H— 

■— i 

ii 

LU 

s: 

Q 

UJ 

cc 

Co 

< 

LU 

s: 

z: 

u 

LL 

o 

c 

o 

• 

o 

X 

CNJ 

o 

II 

*^ 

LU 

o 

X 

LU 
LU  Cl 


(M  CNJ  CNJ  CN 


47 


o 
UJ  jcvj 

00  i 
*v  I 

o 

Q  i 
O 

in 


po  : 
x  o 

o  ; 


O'CT'vOOHof^roLnof^oior^-r- 

r 

j  i 

I 

I  : 

I  i 

i 

orsi(\i^Lr\r\jr\jif\orooooor^QOf\j 

! 


X 

•  i 

to  I 

CD 
< 


r--o>0Lr\irtmir\mm  ir\  irv  ;r\  ir\  un  in 


to 
R 


UJ 

I-  o 
<  LTi 


a  p 


uj  la. 
UJ 

UJ  jCL 


o 

z 

o 

< 

> 

— 

O 

0 

O 

CM 

0 

CM 

LO 

0 

OJ 

fM 

r- 

sO 

sO 

•— 1 

LU 

1— t 

m 

fM 

in 

Cr 

CO 

fM 

m 

>t 

CO 

Cr 

O 

00 

cn 

>fr 

m 

m 

m 

in 

LO 

LO 

LO 

>* 

m 

CM 

•V 

? 

2T 

—J 

(  \ 

V— ' 

X 

»«-*• 

•> 

0 

UJ 

Z 

O 

2: 

a: 

1— 1 

0 

O 

0 

>t 

0 

0 

sO 

>fr 

m 

u 

3E 

r- 

CO 

r- 

O 

<? 

O 

O 

CM 

cr 

>t 

—* 

CVJ 

<M 

CM 

CM 

fM 

CM 

r\j 

fM 

• 

< 

0 

UJ 

a 

X 

«■» 

II 

UJ 

ro 

! 

t— 

2: 

O 

CO 

m 

*0 

O 

CO 

r» 

en 

O 

CM 

ro 

OJ 

i_ 

O 

CO 

r- 

CO 

r- 

LO 

CC 

cr 

cr 

X 

UL 

2. 

• 

0 

• 

• 

• 

• 

• 

« 

• 

• 

• 

0 

O 

O 

CO 

cr 

CO 

r- 

r» 

t—l 

O 

— ' 

• 

LU 

00 

h- 

21 

X 

H- 

X 

> 

O 

D 

0 

a 

_J 

O 

0 

a. 

Z 

O 

O 

oc 

a 

- 

O 

cr 

O 

«=! 

CO 

cm 

fM 

CM 

r- 

r— I 

LO 

r»- 

0 

a: 

h- 

•— 1 

CT 

CM 

CO 

cr 

vD 

m 

CO 

r- 

CC' 

CO 

< 

O 

_J 

• 

• 

• 

• 

• 

a 

• 

9 

• 

• 

• 

• 

a 

UJ 

O 

X 

CO 

CT 

CO 

r- 

r» 

z 

X 

< 

1— t 

00 

• 

-J 

0 

00 

1 

u 

CD 

O 

< 

O 

0 

-J 

O 

0 

lo 

CM 

r- 

0 

r~ 

0 

cc 

fH 

r- 

11 

• 

LU 

O 

m 

» 

CO 

CO 

CO 

X 

r-4 

oi 

2: 

t- 

on 

0 

t— » 

■ — ■ 

1— • 

1—1 

m 

>J- 

m 

m 

m 

CM 

CM 

1— 1 

f— 1 

Q£ 

00 

O 

z 

11 

U_ 

LU 

! 

0 

Z: 

Q 

< 

z 

UJ 

LU 

< 

oc 

ce: 

< 

-J 

ZD 

3 

Q 

a 

f— 1 

1— 

O 

in 

fM 

cr 

CO 

r- 

—* 

fM 

CO 

O 

oo 

2: 

< 

•-4 

r— 1 

0 

cr 

1—4 

a 

r*- 

cr 

fM 

r- 

< 

LU 

< 

i>k 

LU 

_J 

00 

0 

LU 

t— 1 

fM 

LO 

m 

rn 

CM 

— 1 

p— ( 

0 

2. 

«->l 

«— 

If— 1 

Q. 

1 

0! 

UL 

0 

5" 

LL 

a 

LU 

;  1 

LU 

C 

— t 

CX 

X 

• ' 

O- 

00 

LU 

t— 

°1 

oc 

X 

1 

O 

UJ 

;< 

O 

0 

'»! 

LU 

0 

O 

1—4 

£ 

< 

i 

Z 

■=s3 

LU 

0! 

3: 

LU 

Q£ 

CO 

cr 

O 

fM 

m 

LO 

*0 

cr 

Cr 

X 

rsl  1 

i 

LU 

l-H 

.— 1 

l_ 

LU 

O. 

2- 

I 

5: 

1/) 

•— » 

LU 

1 

| 

»— 

O 

4M 


APPENDIX  D-AVERAGE  HOURLY  EVAPORATION 


This  appendix  reports  average  hourly  evapotrans- 
piration  by  the  various  meteorological  equations 
using  the  smoothed  profile  data  and  average  hourly 
radiation  data.  Also  presented  are  the  computed 
Bowen  Ratio  (B)  and  Richardson's  Number  (Rl). 


The  values  in  the  probable  error  columns  gave  the 
range  of  probable  error  that  can  be  expected  in  each 
computed  value  on  the  basis  of  instrument  accuracy 
and  sensitivity.  Missing  data  are  due  to  equipment 
malfunction,  as  indicated  by  an  asterisk. 
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